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Abstract 


An analysis is presented of the effects of airfoil thickness and mean aerodynamic 
loading on boundary-layer receptivity in the leading-edge region. The case of acoustic 
free-stream disturbances, incident on a thin cambered airfoil with a parabolic leading 
edge in a low Mach number flow, is considered. An asymptotic analysis based on large 
Reynolds number is developed, supplemented by numerical results. The airfoil thickness 
distribution enters the theory through a Strouhal number based on the nose radius of 
the airfoil, S = u>r n /U, where u is the frequency of the acoustic wave and U is the mean 
flow speed. The influence of mean aerodynamic loading enters through an effective angle- 
of-attack parameter n, related to flow around the leading edge from the lower surface 
to the upper. The variation of the receptivity level is analyzed as a function of S', fj., 
and characteristics of the free-stream acoustic wave. For an unloaded leading edge, a 
finite nose radius dramatically reduces the receptivity level compared to that for a flat 
plate, the amplitude of the instability waves in the boundary layer being decreased by 
an order of magnitude when S = 0.3. Modest levels of aerodynamic loading are found 
to further decrease the receptivity level for the upper surface of the airfoil, while an 
increase in receptivity level occurs for the lower surface. For larger angles of attack 
close to the critical angle for boundary layer separation, a local rise in the receptivity 
level occurs for the upper surface, while for the lower surface the receptivity decreases. 
The effects of aerodynamic loading are more pronounced at larger values of S. Oblique 
acoustic waves produce much higher receptivity levels than acoustic waves propagating 
downstream parallel to the airfoil chord. 



Chapter 1 
Introduction 


The accurate prediction of boundary-layer transition and the need to control this phe- 
nomenon are becoming increasingly important in the development of advanced aeronau- 
tical vehicles. These requirements have focused renewed attention on the physics of the 
transition process. It was known for many years that boundary-layer transition is influ- 
enced significantly by the free-stream disturbance environment. However, the physical 
mechanisms by which energy is transferred from the long wavelength free-stream distur- 
bances to the much shorter wavelength instability waves were not understood. This came 
to be known as the receptivity problem (Morkovin 1969). 

In order for an external disturbance to generate an instability wave, energy must be 
transferred to the unsteady motion in the boundary layer at an appropriate combination 
of frequency and wavelength. To simplify the discussion, it is useful to consider the 
situation where the external disturbance is of small enough amplitude that the unsteady 
motion can be represented as a linear perturbation of the mean flow. Attention can then 
be restricted to a single time harmonic, with the results for general time dependence 
obtained by superposition. In this report, we consider linear, time-harmonic disturbances 
to the mean flow. 

The earliest theoretical analysis of instability wave generation in a boundary layer 
was presented by Gaster (1965). He considered the case of two-dimensional Tollmien- 
Schlichting (TS) wave excitation in a parallel boundary layer by a time-harmonic dis- 
turbance at the wall. The wall disturbance was localized in the streamwise direction, 
and hence the wavenumber spectrum of this disturbance was broad. Thus, the input 
disturbance contained energy at the appropriate frequency-wavelength combination to 
directly excite an instability wave. 

The first attempts to predict receptivity to naturally occurring free-stream distur- 
bances (Rogler & Reshotko 1975, Tam 1981) were based on a parallel flow formulation 
similar to that of Gaster. However, naturally occurring free-stream disturbances (sound 
waves, turbulence, etc.) travel at much higher speeds than instability waves. Thus, the 
wavenumber spectrum of the free-stream disturbance at a given temporal frequency is 
concentrated at wavenumbers which are significantly different than the wavenumber of 
the instability wave. Hence, these parallel flow analyses succeeded only in finding ‘par- 
ticular solutions’ in the boundary layer which are unrelated to the instability wave. In 
order to transfer energy from a naturally occurring free-stream disturbance to an insta- 
bility wave, a wavelength conversion process is required. A more complete discussion 
of the differences between instability wave generation by localized disturbances and by 
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naturally occurring disturbances can be found in Kerschen (1989). 

Experimental investigations in the 1970s showed that receptivity can occur in the 
vicinity of the leading edge and in localized downstream regions. Reviews of these and 
more recent experiments are presented by Kachanov, Kozlov & Levchenko (1982) and 
Nishioka & Morkovin (1986). The experimental results stimulated theoretical investiga- 
tions (Goldstein 1983, 1985) which showed that the wavelength conversion process takes 
place in regions of the boundary layer where the mean flow exhibits rapid changes in the 
streamwise direction. This occurs (a) near the body leading edge and (b) in any region 
farther downstream where some local feature forces the boundary layer to adjust on a 
short streamwise length scale. The rapid streamwise adjustment requires that nonparallel 
mean flow effects be included at leading order, in contrast to the parallel flow assump- 
tion of classical Orr-Sommerfeld stability theory. Reviews of receptivity are presented by 
Goldstein & Hultgren (1989) and Kerschen (1990). A discussion of the role of receptivity 
in transition prediction can be found in Heinrich, Choudhari & Kerschen (1988). Results 
for receptivity in localized downstream regions are presented in Kerschen, Choudhari & 
Heinrich (1989), Choudhari & Kerschen (1989, 1990) and Choudhari (1990), in addition 
to publications referenced above. 

In this report, we present a theoretical analysis of receptivity in the leading-edge 
region. Previous theoretical studies of leading-edge receptivity have been restricted to 
the Blasius boundary layer on a semi-infinite, zero-thickness plate. The motion has 
been assumed incompressible and two-dimensional. Goldstein (1983) developed a high 
Reynolds number asymptotic analysis for this problem. The small parameter e utilized 
in his analysis is defined as e = (lju/U 2 )~o , where u is the frequency of the unsteady 
motion, v is the kinematic viscosity of the fluid and U is the speed of the free stream. 
Note that e 6 is the frequency parameter F of classical stability theory, which can also 
be interpreted as the inverse of a Reynolds number based on the convective wavelength 
U/u of the unsteady motion. 

Goldstein, Sockol & Sanz (1983) calculated the receptivity coefficient for an acous- 
tic wave propagating downstream parallel to the surface of the flat plate. Heinrich & 
Kerschen (1989) calculated receptivity coefficients for a wide range of free-stream distur- 
bances including obliquely incident acoustic waves, convected gusts (the linear represen- 
tation of free-stream turbulence) of various orientations, and a von Karman vortex street 
passing above the plate surface. At low Mach numbers, the strong leading-edge diffrac- 
tion field created by obliquely incident acoustic waves produces receptivity coefficients 
which are an order of magnitude larger than those for the other types of disturbances. 
The receptivity coefficient for oblique acoustic waves is proportional to M~ a in the low 
Mach number limit, and the strongest receptivity occurs for acoustic waves propagating 
upstream toward the leading edge along one surface of the plate. Heinrich & Kerschen 
(1989) also examined the influence of wind tunnel walls on leading-edge receptivity to 
upstream propagating acoustic waves. A weakly damped oscillatory pattern as a function 
of frequency was found. The oscillations are related to the alternate cut-on of upstream 
and downstream traveling acoustic modes in the wind tunnel. For high frequencies, 
this pattern gradually approaches the isolated plate result. Various parts of Heinrich & 
Kerschen’s results on leading-edge receptivity are summarized in publications referenced 
above; a full account can be found in Heinrich (1989). 

The theoretical investigations of leading-edge receptivity discussed above have signif- 
icantly advanced our understanding of the receptivity process. However, the results are 
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restricted to the incompressible, zero-pressure-gradient Blasius boundary layer on a plate 
of zero thickness. Bodies of practical interest generally have a leading edge with a finite 
radius of curvature, r„, which enters the receptivity problem through the non-dimensional 
parameter 5 = ur n /U. This parameter is typically 0(1) or larger in aerodynamic ap- 
plications. Even in the case of laboratory experiments, 5 may be large enough to cause 
significant deviations from the zero-thickness model. 

A second disadvantage of the zero-thickness geometry is that the influence of mean 
aerodynamic loading cannot be investigated, since any asymmetrical mean flow com- 
ponent in the vicinity of the leading edge leads to boundary-layer separation. It is well 
known that aerodynamic loading significantly influences both the mean development and 
the stability of the boundary layer. Hence, it is reasonable to anticipate that aerody- 
namic loading is an important parameter in the receptivity problem as well. In fact, in 
his laboratory experiments, P. Leehey found (private communication) that small changes 
in the ‘angle of attack’ of his plate produced dramatic changes in the receptivity to acous- 
tic waves. Thus, in order to make predictions of relevance to practical applications, it is 
essential to incorporate both finite leading-edge thickness and mean aerodynamic loading 
into the receptivity theory. 

In the present report, we extend the theory for leading-edge receptivity to address 
the effects of leading-edge thickness and mean aerodynamic loading. In chapter 2, a high 
Reynolds number asymptotic analysis (e <§C 1) of leading-edge receptivity is presented 
for the case of a symmetric mean flow past a parabolic leading edge. The nose radius 
parameter, 5 = ur n /U , is assumed 0(1). We find that, when 5 reaches the value 0.3, the 
receptivity coefficient is reduced by an order of magnitude compared to the 5 = 0 case. 
The numerical results calculated under the assumption that 5 = 0(1) suggest that the 
receptivity coefficient is singular for small values of 5. In chapter 3, the small-5 behavior 
is investigated further by developing a small-5 asymptotic theory. The leading term cor- 
responding to 5 = 0 is found to agree with the flat-plate results of Heinrich and Kerschen 
(1989). For free-stream disturbances which produce a locally symmetric unsteady motion 
about the leading edge, the receptivity coefficient varies linearly for small 5. However, 
for free-stream disturbances which induce an anti-symmetric component in the unsteady 
flow about the leading edge, the receptivity coefficient varies with 52 . This singular term 
leads to a sharp decrease in the amplitude of the receptivity coefficient relative to the 
flat-plate value. In chapter 4, the effects of mean aerodynamic loading are analyzed. The 
influence of aerodynamic loading enters through an effective angle-of-attack parameter 
/i, related to flow around the leading edge from the lower surface to the upper. Modest 
levels of aerodynamic loading are found to decrease the receptivity level for the upper 
surface of the airfoil, while the receptivity is increased for the lower surface. For larger 
angles of attack close to the critical angle for boundary layer separation, a local rise in 
the receptivity occurs for the upper surface, while on the lower surface the receptivity 
decreases. These effects are more pronounced at larger values of 5. Conclusions of the 
study and recommendations for future work are summarized in chapter 5. This work has 
also been presented in Hammerton & Kerschen (1996, 1997, 1998). 
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Chapter 2 


Receptivity for a leading edge with 
no aerodynamic loading 

2.1 Introduction 

The receptivity process through which free-stream disturbances generate instability waves 
in boundary layers was first discussed by Morkovin (1969). The transfer of energy from 
the free-stream disturbance to the instability wave generally comes about through non- 
parallel mean flow effects, which may arise either in the leading-edge region, or in a 
localized region farther downstream in the boundary layer (Goldstein & Hultgren 1989; 
Kerschen 1990). 

Up to now, theoretical studies of leading-edge receptivity have been restricted to 
a semi-infinite, zero-thickness plate. Goldstein (1983) developed an asymptotic analysis 
for this problem; leading-edge receptivity coefficients for various free-stream disturbances 
were calculated by Goldstein et al. (1983) and Heinrich & Kerschen (1989). However, 
aerodynamic bodies designed for subsonic flow generally have finite thickness distribu- 
tions with a parabolic leading edge. In the present chapter, we examine the influence 
of the thickness of a body on leading-edge receptivity. The body is assumed to be two- 
dimensional, with a symmetric cross-section and a parabolic leading edge. An asymptotic 
theory for the case of a symmetric mean flow is developed in this chapter. Results are 
presented for receptivity to acoustic waves in the free stream, incident on the body at 
arbitrary angle. Attention is focused on the variation of the receptivity level with the 
nose radius of the body and the incidence angle of the acoustic field. 

In §2.2, a high Reynolds number asymptotic analysis (e 6 = uu>/U 2 \ e <g: 1) is formu- 
lated for an incompressible, two-dimensional flow. This follows the approach of Goldstein 
(1983) but remains valid for a nose radius comparable to the free-stream disturbance 
length scale U e /u>. Here U e is an effective free-stream speed, defined in §2.2, which in- 
cludes a correction to the free-stream velocity due to the flow perturbation created by 
downstream portions of the airfoil. Two streamwise regions enter the analysis, one re- 
gion where the distance downstream is 0(U e /ui) and the disturbance is governed by the 
linearized unsteady boundary layer equation (LUBLE), and a second region at distances 
0(e~ 2 U e /u) where the disturbance is governed by the triple-deck structure, correspond- 
ing to the asymptotic form of the Orr-Sommerfeld equation (OSE) in the vicinity of the 
lower branch. 

The structure of the solution in the LUBLE region is analyzed in §2.3. The inviscid 
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pressure field and slip velocity induced by the free-stream disturbance (analyzed in §2.4) 
drives the unsteady motion in the boundary layer. Far downstream in the LUBLE region, 
the solution consists of a Stokes wave, and a set of asymptotic eigensolutions that contain 
velocity but not pressure fluctuations. These asymptotic eigensolutions are equivalent to 
those obtained by Lam & Rott (1960) and Ackerberg &: Phillips (1972), but modified to 
include the effects of the mean pressure gradient and surface curvature. The asymptotic 
analysis for distances far downstream determines the form of the eigenfunctions, but 
not their coefficients C*. These coefficients can be found only through a full solution 
of the LUBLE, which must be determined by numerical methods. Lam & Rott (1993) 
have recently generalized their eigenfunctions to take account of arbitrary streamwise 
variation in the mean flow. The direct development in parabolic coordinates presented 
here is more convenient for our purposes, but it can be shown that our expressions for 
the eigenfunctions are in agreement with these more general results. 

The wavelengths of the asymptotic eigenfunctions shorten progressively with distance 
downstream. Eventually, the self-induced pressure field associated with the displacement 
thickness of each asymptotic eigenfunction becomes significant, and the triple-deck struc- 
ture replaces the LUBLE as the correct asymptotic approximation to the Navier-Stokes 
equation. The first asymptotic eigenfunction of the LUBLE matches on to the Tollmien- 
Schlichting wave solution of this triple-deck region. Thus, the form of the free-stream 
disturbance and the geometry close to the nose influence the amplitude of the Tollmien- 
Schlichting wave only through the coefficient C\ of the first asymptotic eigenfunction. 
Therefore, we call C\ the ‘Receptivity Coefficient’. The primary objective of this chapter 
is to determine the Receptivity Coefficient as a function of leading-edge geometry and 
free-stream disturbance characteristics. 

For the flat-plate case, it is not clear whether the set of eigensolutions obtained by 
Lam & Rott are complete. A second, very different set of eigensolutions was obtained by 
Brown & Stewartson (1973) which, they argue, better represent the physical properties of 
the flow. While the relationship between these two sets is a fundamental question which 
deserves further study, this must first be done in the context of the flat-plate problem. For 
a parabolic body, generalizations of the Brown & Stewartson eigensolutions should also 
exist, but we concentrate solely on the generalizations of the Lam- Rott solutions since it 
is demonstrated in Appendix B that they match naturally to the Orr-Sommerfeld modes 
further downstream. 

In §2.4, the inviscid pressure field and slip velocity produced by the interaction of a 
free-stream acoustic wave with an airfoil are determined. The exact form of the unsteady 
slip velocity in the vicinity of the leading edge is determined by the global solution about 
the airfoil. This depends on the magnitude of the reduced acoustic frequency k = ub/c, 
where b is the airfoil semi-chord and c is the speed of sound. Here we present results 
for the limiting cases k < 1 and k » 1, when relatively simple expressions for the slip 
velocity can be obtained. Numerical solutions of the LUBLE are then carried out in §2.5, 
and comparisons with the asymptotic eigenfunctions of §2.3 are utilized to determine the 
Receptivity Coefficient C x as a function of S = ur n /U c and characteristics of the free- 
stream acoustic wave. The analysis presented is for S = 0(1); the only restriction is that 
S < e -2 , so that the unsteady disturbance in the nose region is governed by the LUBLE. 

2.2 Formulation 
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Figure 2.1: An illustration of the physical situation of interest: a thin, symmetric airfoil 
of chord 2 6 is at zero angle-of-attack in a uniform flow of speed U Q 0 , with a plane wave 
incident at an angle 0 with respect to the airfoil chord. 


We consider a thin, symmetric airfoil of chord 26 at zero angle-of-attack in a uniform 
flow of speed C/^,. A plane acoustic wave of frequency u > , propagating at an angle 6 
with respect to the airfoil chord, is assumed to be incident on the airfoil as illustrated 
in figure 2.1. Two-dimensional, low-Mach-number flow is considered. Since the Mach 
number is small, the mean flow can be analyzed using incompressible theory. For the 
unsteady component of the flow, most features of interest can also be analyzed with 
incompressible theory. The influence of compressibility on the unsteady component of 
the flow is discussed in §2.4. The Reynolds number is assumed large, so the flow field is 
inviscid and irrotational everywhere except in the vicinity of the airfoil surface. 

2.2.1 Inviscid Outer Flow 

Introducing Cartesian coordinates (x, y ) normalized by the airfoil semi-chord 6, with the 
origin located at the airfoil leading edge, the airfoil shape is given by 

y = ±6s(x), 0 < x < 2 (2.1) 

where s (x) is the non-dimensional thickness distribution of the airfoil and the thickness 
parameter <5 « 1. The airfoil is assumed to have a rounded leading edge and a sharp 
trailing edge. The dimensional complex potential W for the steady, inviscid flow past 
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(2.2) 


the airfoil is given by thin-airfoil theory, 

W = U 00 b z + 8— f 2 In (z — X[) da?! + O (5 2 ) 

7T Jo ' ' 

where z = x 4- \y and the dimensional velocity (u, v) is given by u — iu = b~^dW/dz. 

For an airfoil with a rounded leading edge of radius r n , the small argument expansion 
of the thickness distribution has the form 

s ( x ) = SjX I/2 + s 2 x + s 3 : r 3/2 + 0(x 2 ) (2.3) 

where = (2 r n /8 2 b) x ^ 2 . The coefficient of the leading term in 2.3 must be O (1), implying 
that r n = 0(8 2 b). The thin-airfoil expansion 2.2 is invalid in the vicinity of the leading 
edge, and must be replaced by a local expansion (Van Dyke, 1964a: Chapter 4). The 
natural length scale for this local expansion is the airfoil nose radius r n . Introducing a 
parabolic coordinate system ffj in the local region, with r n as the length scale, 

x+iy= \t [(^ +i ^) 2+1 l ’ ( 2 - 4 ) 

the airfoil surface is given by fj = 1 + |<fs 2 ^ 3 /(l + £ 2 ) + 0(8 2 ) . Thus, at leading order in 
the local coordinates, the airfoil surface is defined by the parabola fj = 1. 

The complex potential describing the leading approximation to the steady, inviscid 
flow in the vicinity of the airfoil nose is 

W = \u c r n (£+* (fj- l)) 2 , (2.5) 

leading to a slip velocity, f/ e £/(£ 2 + l) 1 ^ 2 , where the constant U e is determined by match- 
ing with the thin airfoil expansion (2.2). For leading-edge shapes that contain a wedge 
component (S 2 ^ 0), the approximation (2.5) is valid only at 0(1) in the thickness pa- 
rameter 8 and matching with (2.2) then shows that U e = U^. However, for leading-edge 
shapes in which the S 2 multiplies a locally analytic function of x, the even coefficients 
in (2.3) vanish. The Joukowski airfoil is one example of such an airfoil. In this case the 
leading-edge region is also parabolic at 0(5), and the matching then gives 



The 0(8) term in (2.6) is essentially a correction to the ‘free-stream speed’ in the local 
leading-edge region, due to the flow perturbation created by downstream portions of the 
airfoil. 

The scattering of the acoustic wave by the airfoil produces an unsteady perturbation 
to the inviscid, irrotational flow described above. This unsteady, inviscid perturbation is 
considered in §2.4. The slip velocity and pressure associated with the inviscid flow drive 
the viscous flow in the boundary layer adjacent to the airfoil surface. 

2.2.2 Boundary-layer flow 
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The boundary-layer flow in the vicinity of the leading edge is also analyzed most conve- 
niently in parabolic coordinates. Although the nose radius r„ is the most natural length 
scale for the steady flow, the length scale U c /u is more convenient for analysis of the 
unsteady flow. Thus, we introduce new coordinates, 


T} = S?fj, 


(2.7) 


where 



(2.8) 


is a Strouhal number based on the airfoil nose radius. 

The flow around the body is analyzed in terms of the incompressible vorticity equa- 
tion. Utilizing U e for the velocity scale, U e /u for the length scale and u~ l for the time 
scale, the non-dimensional form of the vorticity equation is 


2 an a(n,'l') 6 /a 2 n a 2 n\ 

ft + ~ e (aja + dij 2 ) 


(2.9) 


where ^ is the streamfunction, the velocity is given by v = h _1 (^^e ? - — ty/e^), where 
e £ and are unit vectors in the f and fj coordinate directions, respectively, and the 

metric coefficient h = (£ 2 4- fj 2 )* . The vorticity —Q is related to the stream function by 
n = h~ 2 (4>^ + tyrjrj)- The small parameter 


c UUJ 

e =yT « 1 ( 2 W) 

is the reciprocal of the Reynolds number based on the disturbance lengthscale. The flow 
satisfies the no-slip boundary conditions on the body surface, 

tj> = tyj = 0, on 77 = 52, (2.11) 

where terms of higher order in S have been neglected. The analysis presented here 
considers the small-e limit but with 5 = 0(1). 

Since the Reynolds number is assumed large, viscosity is important only in a thin 
boundary layer adjacent to the body surface. To analyze the boundary-layer flow, we set 

fj - 52 = e 3 77, | = f = e 3 #. (2.12) 

The incompressible vorticity equation, expressed in terms of the stream function, then 
becomes 






\b 

H 2 


* 




H 2 


+e 


12 


+ 








H 2 






[H 2 \ 




d(^V) 



(2.13) 


where H — (£ 2 + 5)2. This equation is exact, except for the approximation that h = H in 
the boundary layer. Careful consideration of the asymptotic structure, at all stages of the 
evolution of the disturbance, shows that the correction term 0(e 3 r]H~ 2 ) never becomes 
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significant and hence it is dropped throughout this presentation. At large values of 77, 
the boundary-layer flow matches to the inviscid slip velocity, 

— 1 ‘ U a (€,t), as 77 — *■ 00. ( 2 . 14 ) 

For 0(1) values of £, the terms on the right side of (2.13) can be neglected, leading 
to the unsteady boundary-layer equation. The unsteady component of the flow, a small 
perturbation to the mean flow, then satisfies the linearized unsteady boundary layer 
equation (LUBLE). The LUBLE region is considered in the following section. However, 
the solution of the LUBLE contains components whose wavelengths progressively shorten 
with distance downstream. When f = 0(e _1 ), terms on the right hand side of (2.13) 
become significant and the correct asymptotic approximation to (2.13) has the triple- 
deck structure. The structure of the developing boundary layer is summarized in figure 
2.2. The asymptotic matching of the two streamwise regions is discussed briefly at the 
end of §2.3. 

2.3 Unsteady boundary- layer region 

In this section, we consider (2. 13) for 0( 1 ) values of f in the limit e -► 0. Subsequently, the 
behavior of this solution at large values of £ will be considered, in anticipation of matching 
with the Orr-Sommerfeld region that exists farther downstream in the boundary layer. 
F° r £ = 0(1) and e — ► 0, the terms on the right-hand side of (2.13) can be neglected. 
Integrating once with respect to 77 then gives the unsteady boundary-layer equation 

H 2 V vt + (%#€„ - *77^) - ZH-Hf -4V, = -H 2 Ps (2.15) 

where the pressure gradient P ? (f,t) is obtained by matching to the outer (inviscid) flow, 

In this chapter, we consider uncambered airfoils at zero angle-of-attack to the mean 
flow. The slip velocity in the region of the parabolic nose is then given by 

U 3 (S,t) = jj + Au a (Z)e- u , (2.17) 

where the steady contribution follows from (2.5), while the time-dependent perturbation 
depends on the particular form of the free-stream disturbance, as discussed in §2.4. Since 
the velocity field associated with an acoustic wave is of very small amplitude compared 
to the mean-flow speed, we assume that the amplitude scale factor A « 1. The steady 
and unsteady components of the flow field can then be analyzed separately. Thus, the 
streamfunction within the boundary layer can be written in the corresponding form 

V = £0(£,q) + A^(£,77)e -i '. (2.18) 

where extraction of the factor f from the mean-flow component is motivated by the form 
of (2.17). The function 0(£, 77) describing the steady boundary-layer flow satisfies 

^7 vv 4" ^77^ T £(0770^ — 0^0^) SH (0 fJ 1 ) = 0 , (2.19) 
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u* 


Figure 2.2: A schematic illustration of the boundary layer structure for a body with a 
parabolic leading-edge in a symmetric mean flow. The three decks in the Orr-Sommerfeld 
region are (i) the viscous wall layer; (ii) the main inviscid layer; and (iii) the outer 
irrotational layer. 
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with boundary conditions <£ = <&, = 0 at 77 = 0 and 4> n — * 1 as 77 — > 00. It may be 
observed that written in this form the steady equation involves S and hence it appears 
that the steady flow is dependent on u>, the frequency of the unsteady perturbation. 
However, u> also enters the scaling of the streamwise coordinate £. Re-writing (2.19) in 
terms of £ = S~^, 

<t>rm + K<i> + - <MV-) - (1 + f )" W - 1) = 0, (2.20) 

it is seen that the steady flow is indeed independent of uj. 

The time-dependent contribution to the boundary-layer flow satisfies the linearized 
unsteady boundary layer equation (LUBLE), 


T 2 d P 


HV) = VWl +(<£ + + ( i//2 - + ~ H 2 S <t>r, 

du 3 
H d£ ’ 




with boundary conditions 

ip = ip n = 0 at 7 ? = 0 and xp v — ► Hu s 


( 2 . 21 ) 

as 77 — ► 00. (2.22) 


These equations must be solved numerically; the methods used are described in §2.5. 
However, as discussed in §2.1, we are primarily interested in the component of the un- 
steady field that, in the large-£ limit, matches onto the Tollmien-Schlichting wave so- 
lution of the Orr-Sommerfeld equation. In the following two subsections, we develop 
large-£ asymptotic solutions for the steady and unsteady components of the flow. For 
convenience, in the rest of this section we restrict attention to the upper surface of the 
body (ie £ > 0). Corresponding results for £ < 0 can be obtained by inspection. 


2.3.1 Steady Boundary Layer Equation 

In order to analyze the large-£ behavior of the unsteady flow, the asymptotic form of the 
steady flow is required. The expansion for the steady flow is most naturally developed in 
terms of the streamwise coordinate scaled on the body nose radius, £, as defined in (2.4). 
Far downstream of the nose (£ ;» 1), the pressure gradient for the steady flow past a 
parabola decays to zero and hence the mean flow approaches the flat-plate solution. The 
asymptotic form of the solution is given by Van Dyke (1964b), 

0(f . n) ~ F(, 7 ) + + (B.C.ft) + C,(i/))i, (2.23) 

where F( 77) is the Blasius function, Gi (77) = tjF'(t 7) - F(r/), A\ and B\ are as yet unde- 
termined numerical constants, and £2(77) satisfies 

G 2 + F G 2 +2F'G 2 — F"G 2 = F ,2 -l+2A l FF" ) G 2 ( 0) = G' 2 ( 0) = C"( 0) = 0. (2.24) 
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The value of A\ is fixed by the physical requirement of exponential decay of vorticity at 
the outer edge of the boundary layer. Hence G 2 {t}) must decay exponentially for large 77, 
which gives A\ = 0.60115. The value of B\ cannot be determined by the large-£ analysis 
since £,~ 2 G\(rj) is an eigensolution of the perturbation equation. The next term in the 
expansion is 0(£ _Q ), a « 3.774, the fractional power arising as the next eigensolution 
of an infinite sequence (Libby & Fox, 1963). It appears that B u together with the set 
of similar constants appearing in higher-order terms, is dependent on conditions close to 
the nose of the body and hence can be determined only by numerical integration from 
£ = 0. This we discuss in §2.5. The next four ignored terms in the expansion are of order 
£-3.774^ £ -4 ln<f and f -4 . As we will demonstrate in due course, the fact that 

these terms are of similar magnitude until <f is extremely large poses certain problems in 
the numerical treatment of the problem. 

In developing the solution to the LUBLE, the limiting behavior of the steady flow 
close to the surface is required. This is found to be 






2 


^1 4 - A\ 


In £ 2 

T~ 



rf 1 
~6? 


Up V 

5 ! 


+0tfr 3 - 77 Wr 2 Ino, 


(2.25) 


where Uq = F"(0) = 0.4696. 


2.3.2 Linearized Unsteady Boundary Layer Equation 

We now consider the evolution of the unsteady perturbation to the mean flow, which is 
governed by (2.21). Far downstream (f > 1), the unsteady component of the stream- 
function consists of a particular solution, determined entirely by the local conditions 
far downstream, together with a set of asymptotic eigensolutions, 

6 S) = Mv, £;$) + £ Gi(S)ipi(r], e S). (2.26) 

t 

The particular solution is a generalization of the classical Stokes layer solution, driven by 
the local value of the unsteady pressure gradient (2.21), see Lighthill (1954) for details. 
The eigensolutions ipi depend on the geometry of the body far downstream, but are 
independent of the local free-stream disturbance. These eigensolutions are generalized 
forms of the eigenfunctions found by Lam & Rott (1960) for the flat-plate boundary 
layer, taking account of the non-Blasius mean flow. The coefficients C* multiplying the 
asymptotic eigenfunctions are determined entirely by conditions close to the leading edge 
(f = 0(1)). One of these eigensolutions, which we label tpi, matches on to the Tollmien- 
Schlichting wave in the Orr-Sommerfeld region farther downstream, where £ = 0(e~ l ). 
Thus, it is only through the coefficient C\ that the unsteady disturbances in the free 
stream influence the amplitude of the Tollmien-Schlichting wave. 

Our primary interest is in the relationship between the free-stream disturbances and 
the amplitude of the Tollmien-Schlichting wave. Thus, we focus on the asymptotic eigen- 
solutions of the LUBLE. Guided by the results for the flat-plate case (Goldstein, 1983), 
we anticipate the development of a two-layer structure for f ;» 1, consisting of a main 
layer where 77 = 0(1) and a new inner layer of width 0(£ -1 ). 

Inner layer 

Introducing an inner variable m = £77, the homogeneous form of the LUBLE then becomes 
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(2.27) 


”071 


+ K 1 + 


m + ~ mm ) + 0 



^ j^o 7(C/^)(m^ - - 2^) - - 4^)| 

where it follows from (2.23) that 7(f) = 1+2^! ln£/f 2 + B v /^ + 0{l~ ZTlA ). The regular 
perturbation solution to (2.28), for large £, leads to solutions related to the particular 
solution ipp. Thus, the asymptotic eigensolutions must arise as a balance between the 
highest m derivative ( i.e . wall-normal derivative) and the terms involving £ derivatives. 
It is then clear that the eigensolution must contain a factor exp(T(£)) and it is for this 
reason that the higher order terms containing derivatives with respect to £ have been 
retained in (2.28). 

Writing ip = exp (T(£))/(£, M), where a new boundary layer variable M = m(l + 
S/2^ 2 ) has been introduced to simplify the solution, /(£, M) satisfies 

h,MM + ifM - - | /VS) |(M/« -f) = Ji R + °(p. (2.28) 

where f(£) = 7(f) - 3/2£ 2 and 


ft ~ UtefMfa — ft) + ~^(M 2 /mm — 2MJm) 


UfM 3 d T 
4!£ 2 df 


(M/m - 4/) - 


SMdT 


(M/ m - 2/). 


(2.29) 


If dT/d£ is set equal to — \£ 2 /UqT, where A is an eigenvalue to be determined sub- 
sequently, solution to the leading order equation immediately follows using separation of 
variables, 

if ~ D? T e T M [p(M) + C 3 q(M) + ...]. (2.30) 

Here D is an arbitrary constant whose value is chosen for convenience later, and 7 is a 
constant that is determined at a later stage of the analysis. The exponent T(£) is then 
determined by integrating the large-£ expansion of f 2 /r(f/S2). This exponent can be 
separated into two components, T = T^ 0) + T^ l \ where 


T (0) (O 


xe 

Uh 



ne/s) 

e 


+ (2Ai + 3/2 — B\) 



(2.31) 


is large for large £ and hence must be included in the leading-order form of the eigenso- 
lution, while the terms in TO)(£) are small for f 1 and thus could be considered as 
algebraic correction terms rather than as exponential amplitude functions. The first few 
terms in are of the form 


T<‘>({) = + t 2 5 2 ' n2( ^ /S> + + 6 4 S 2 i + O (r 2 ' 635 ) , (2.32) 
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where the coefficients b, are determined by the higher order terms of (2.25). 

The leading-order mode-shape function p(M) satisfies the differential equation 


C{p) = 0, 


(2.33) 


where the operator C is defined by 

C(p) = p'" + ip' + \{Mp' - p), (2.34) 

with p' = p \[ , etc. The boundary conditions at the wall are 


P(0) = p'(0) = 0, 


(2.35) 


and matching to the main layer requires that p" — > 0 as M — * oo so as to avoid expo- 
nential growth of the inner solution. The differential operator and boundary conditions 
are homogeneous and thus form an eigenvalue problem for A. In fact, the leading-order 
equation is identical to the flat-plate equation, the effect of curvature having been ab- 
sorbed into the boundary- layer variable M. Thus p(M) and the value of the eigenvalue 
A = e~* ,n p ~2 can be obtained immediately from Goldstein (1983) as 


gy t'(M - M)Ai(z)dM 
/o°° Ai(5)dM 


p(M) = 


4 i ?r — i i ~r 


z = e4 m p 2 M — p 


(2.36) 


where p is a solution of Ai'(— p) = 0. An infinite set of such roots (pi > 0) exists, with 
corresponding mode-shape functions Pi(M). 

The value of r for each eigensolution, ip, is then determined by a solvability condition 
on the 0(£~ 3 ) correction to the mode shape, q(M). This function satisfies 


C(q) = 7 l(p), 


(2.37) 


where 


n{p)=2TUHMp'-p) + ^(MY-2Mr/)-^-(Mp'-4p)+^-(Mp'-2p). (2.38) 

The boundary conditions at the wall and the matching condition to the main layer are 
the same as for p(M). Integrating by parts and using the boundary/matching conditions 
on p and q gives the relation 

roo roo 

J o p'"C(q)dM = J o q'"C{p)dM, (2.39) 

the right hand side vanishing by virtue of (2.33). Hence r is determined by the condition 

roo 

/ p'"n(p)dM - 0. (2.40) 

J 0 

This gives r = r^ 0) + St^ x \ where r^ 0) is the fiat-plate value and 

(!) £> /o°° MpHJW ~ 2p)dM 

4 Ufpl /o°° P >, '{Mp' — p)dM • 1 ' j 
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It should be noted that this result does not involve linearization in 5, but is valid for 
S = 0(1). By suitable change of variables, this can be re-cast into the form 


i J 2 — 2Ai 2 (— p) 
4 Wp Ti 

8 p 3 — 27 . 

4ot^y 1 


(2.42) 


where expressions for the integrals J n = /“(z + p)"[Ai(z)] 2 dz are given in Appendix 
A. Similarly, the expression for (Goldstein, 1983) can be expressed in terms of p, 
avoiding the need for numerical evaluation of the integrals, 


( 0 ) _ 1 /_ 3 , = 889 - 16p 3 

4 V 12 Jj 1260 ' 


(2.43) 


The value of 77 for each eigensolution is obtained by substituting pi into (2.42) and (2.43). 
Finally, we set the value of the arbitrary constant in (2.30) to D = 2~ T so that, in the 
limit S — ♦ 0 , our expression for the asymptotic eigenfunction reduces to the result given 
by Goldstein (1983). 


Main layer 

In the main part of the boundary layer, where 77 = 0(1), the unsteady motion is es- 
sentially an inviscid response, driven by the displacement thickness of the inner layer. 
Setting rp = (£ 2 / 2 ) T exp(T(f))< 7 (£, 77 ) and substituting into the homogeneous form of the 
LUBLE, g(£, rj) is found to satisfy 

<P m 9 - A,9„ = ei ^ + 0 (r 3 ). (2.44) 

Solving, and matching to the inner layer then gives 

9 = i (f*(!7) + e ‘"^ Ui + 0(r 2 ) ) ■ (2.45) 


Large-^ matching 

In Appendix B it is shown that it is the first asymptotic eigensolution which matches 
on to the Tollmien-Schlichting wave that becomes unstable farther downstream. In §2.5, 
numerical solutions of the LUBLE are compared to the first asymptotic eigenfunction, 
in order to extract the receptivity coefficient C\. The two most convenient points of 
comparison are the wall shear and the displacement thickness of the boundary layer. 
The wall shear is determined by the solution in the inner layer. However, for comparison 
with numerical solutions, it is most convenient to express the wall shear in terms of a 
derivative with respect to 77 . Rewriting (2.30) in terms of 77 and evaluating at 77 = 0, 
the contribution to the wall shear due to the asymptotic eigensolution ipi is 


^"(77 = 0) 



5_ p Ti(Oc2. 


U' Q p- 


2 eV K 


(i + 0 (r 2 )). 


fo°° Ai(z - pi) dz 


(2.46) 


15 



Taking the first root of Ai'(-p) = 0, namely p { = 1.01879, and evaluating the numerical 
coefficient, the wall shear for the first asymptotic eigensolution is 

ipi(V = 0) = 0.4356(1 + i)e r ‘ >(0 e T (,) («> (i + 0(£ -2 )), (2.47) 

where T[°\z) is given by (2.31), Ti^(£) is a small correction term given by (2.32) and 


t { = -0.69213 -S 1.9878 i. 


(2.48) 


The contribution to the displacement thickness due to an asymptotic eigensolution is 
obtained from the solution in the main layer, 


MV 


oo) ~ 



‘ e™£ 



+ o(c 2 ) 


(2.49) 


where the result for the first asymptotic eigensolution ipi follows by substituting pi and 

Since T oc £ 3 at leading order, the wavelengths of the asymptotic eigenfunctions 
decrease with distance downstream, increasing the importance of longitudinal derivatives. 
The LUBLE (2.15) is obtained from (2.13) by neglecting higher-order derivatives with 
respect to £. Thus, although the asymptotic eigensolutions derived in this section are 
uniformly valid solutions of the LUBLE as £ — ► oo, they are not uniformly valid solutions 
of the full linearized Navier-Stokes equation. 

Terms on the right side of (2.13), which were neglected in the LUBLE, become sig- 
nificant when £ = 0(e -1 ). An irrotational layer outside the mean boundary layer, driven 
by the oscillating displacement thickness of the eigensolutions, must then be considered, 
and the associated pressure gradient appears in the leading-order equations governing 
the inner layer. This coupled viscous-inviscid interaction has the triple-deck structure, 
corresponding to the small-e asymptotic approximation to the Orr-Sommerfeld equation 
in the vicinity of the lower branch. The development of the boundary layer structure 
was summarized in figure 2.2. A complete treatment of the linear development of the 
instability would require an asymptotic solution for the Tollmien-Schlichting wave of the 
Orr-Sommerfeld equation, taking account of the surface curvature and non-zero pressure 
gradient. This Tollmien-Schlichting wave solution could then be matched to the first 
Lam-Rott asymptotic eigensolution of the LUBLE, as was done by Goldstein (1983) for 
the flat-plate case. However, for 0(1) values of 5, with e and <5 of the same order, the 
Orr-Sommerfeld region is influenced by the full airfoil thickness distribution, necessitat- 
ing a general development that does not seem justified in the present context. For flow 
around a semi-infinite parabolic body, we verify in Appendix B that the eigensolutions 
obtained here do indeed match on to the Tollmien-Schlichting wave. 

In this chapter we restrict ourselves to examining only the receptivity process, that 
is, the determination of the coefficient C { of the first eigensolution, which matches to 
the unstable Tollmien-Schlichting mode. The remainder of the chapter is concerned 
with obtaining the Receptivity Coefficient for different free-stream disturbances. This is 
accomplished by comparing numerical solutions with the asymptotic solutions obtained 
above. 
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2.4 Free-stream Disturbances 


In this section we consider the inviscid flow field produced by the interaction of a free- 
stream acoustic wave with the airfoil. The acoustic wave is assumed to be incident on 
the airfoil at an angle 9 with respect to the airfoil chord, as illustrated in figure 2.1. The 
slip velocity and surface pressure fields generated by this interaction drive the unsteady 
motion in the boundary layer, leading to the generation of a Tollmien-Schlichting wave. 

For a low Mach number flow, the acoustic wavelength 2ttc/u is long compared to 
the hydrodynamic length scale Uoo/u). Thus, outside the boundary layer, the unsteady 
flow in the vicinity of the leading edge is incompressible and irrotational. Potential flow 
theory then shows that this local flow has the form 

««(0 = *. (0)jj + (2.50) 

Here 9,/H and 1/H correspond to purely symmetric and anti-symmetric flow about the 
leading edge, respectively. The coefficients k 3 and K a multiplying these eigenfunctions 
are independent of the nose geometry, but depend on the free-stream disturbance, being 
determined by global features of the unsteady flow. Substituting (2.50) into (2.21) gives 

T(ip) = k s £ | \H 2 — — j + K a ^i// 2 + - — ^ . (2.51) 

Hence, writing ip = K 3 ip 3 + K a ip a , it follows that the Receptivity Coefficient for the upper 
surface is given by 

Ci(S) = k,{9) C 3 {S) + K a {6) C a (S), (2.52) 

where C 3 and C a are extracted from the solutions of 

Hl>*) = £ (iH 2 
0 F(4>a) = iH 2 + 


_ 2 S\ 

H 2 ) 

£ 2 ~S 

H 2 


(2.53) 


respectively. The corresponding Receptivity Coefficient for the lower surface is obtained 
by replacing 9 by —9. 

The remainder of this section is concerned with calculating /c a ($) and K a (9), the 
coefficients of symmetric and anti-symmetric flow about the nose. These coefficients 
are found by asymptotic matching of the local solution (2.50) for the unsteady slip 
velocity with an appropriate global solution. The nature of this solution depends upon 
the magnitude of the reduced acoustic frequency k = ub/c , where c is the speed of 
sound in the undisturbed medium. Calculation of the velocity field about a thin wing, 
including compressibility effects, is described in Sedov (1965: Chapter 2, pp 87-107). 
In general no simple expression for the slip velocity can be obtained. Here we consider 
two cases that do lead to relatively simple results. The first case is that of extremely 
low Mach numbers, such that the acoustic wavelength is long not only compared to the 
hydrodynamic length scale, but also compared to the airfoil chord, i.e., k <C 1. 

In this situation the unsteady interaction of the acoustic wave with the airfoil can be 
analyzed using the classical unsteady airfoil theory for incompressible flow. The second 
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case is that of acoustic wavelengths long compared to U 0 0 /u but short compared to the 
airfoil chord, i.e., k » 1. In this case acoustic diffraction theory can be used to analyze 
the interaction. 

First consider the case k <?C 1, where incompressible, unsteady airfoil theory applies 
(see for example, Garrick, 1957). Since the airfoil semi-chord is the relevant length scale 
for the unsteady aerodynamic interaction, the airfoil thickness can be ignored (except in 
the region of the leading-edge, where (2.50) is applicable). Thus, the airfoil reduces to 
a zero-thickness flat plate, and the solution consists of two potential flow components, a 
non-circulatory component, plus a component due to the vorticity shed from the sharp 
trailing edge. The shed vorticity is convected downstream by the mean flow and hence 
takes the form j(x, t) = where a = ub/Ue. is the aerodynamic reduced frequency. 

The non-dimensional complex potential for this unsteady flow is given by 





i7e lt 
2n 



,C - 1/Co. 


dx 0 , 


(2.54) 


where £( 2 ) = 2 — 1 + ^ 2(2 — 2) and Co = C( x o)- Applying the unsteady Kutta condition 
then fixes the strength of the shed vorticity, 


7 = 


27rsin£? 


J 2 00 e iax ° 


te)' d *° 


(2.55) 


For matching with the local solution (2.50), the small-x expansion of the velocity on the 
airfoil surface is required. We find 


u(x, 0±) ~ |±^^(1 -f J(a)) x 2 + cos# + 0(x2)j e 14 as x — ► 0, (2.56) 

where J(a) represents the effect of the shed vorticity, 


J(a) = 


(«=»)* dst c m 
jr *“»(*)* d*„ - a ' 


as a — > co. 


(2.57) 


Thus for the high reduced aerodynamic frequencies of interest in the present study, the 
effect of shed vorticity on the unsteady flow in the vicinity of the leading edge is smaller 
than the contribution due to the non-circulatory component, by a factor of 0(l/a). 
Matching with the large £ limit of the local leading-edge solution (2.50) then gives 


k 3 = cos 9 , 


= a2 sin 6 . 


(2.58) 


The symmetric and anti-symmetric components of the local flow past the leading edge 
are seen to be in phase, owing to the incompressible nature of the unsteady interaction. 
Except for values of 9 near 0 and 7r, the slip velocity near the leading edge is dominated 
by the anti-symmetric component of (2.50), which reaches a peak value (b/r n ) 1 / 2 sin 9 at 
the nose of the airfoil. Note also that, in this limit, the flow about the nose (and hence 
the Receptivity Coefficient) is the same for an airfoil with sharp trailing edge as for a 
body with a rounded trailing edge. 
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Next consider the opposite limiting case, k » 1. The acoustic wavelength is then 
short compared to the airfoil chord, and the interaction of the acoustic wave with the 
leading edge can be analyzed by taking the airfoil chord to be semi-infinite. As in the 
first case, on the scale of the acoustic wavelength, the airfoil appears at leading order as 
a zero-thickness plate and the problem reduces to the classical Sommerfeld diffraction 
problem (Noble, 1958). The velocity potential for the unsteady flow is given by 


<t>{x,y) = |-i e ifc(lcose+J/sine) - sin ±9 
which leads to a slip velocity 


• i a sgn(y) [°° exp(-(A 2 - A: 2 ) i |y| - iAx) 

~ ~in f / 

k 7 r( e lk\2 J - 


2 v 1 i t — - dA > e , 

7i(2k) 2 J-oo (A + k cos 6) ( A + k) 2 

' (2.59) 


a 1 


u(x, 0±) = l cos9e' kcosd (l erf$) ± sin \9 


where 4> = e *' n (2kx) 2 sin ±9. Hence, as x — ► 0, 




-u 


(nk) 


x 2 > e 


2 2 e ~ 4 57r I 

u(x, 0±) ~ l ± sin ~ 6 1“^”^ + cos 9 + 0(x 2 ) > e _lt . 

(7rA:) 2 


(2.60) 


(2,61) 


Matching the slip velocity on the upper surface to the local solution in the vicinity of the 
leading edge, then gives 


K s = cos 9, 


Ka 


2e~« i,r 


sin {\9), 


(2.62) 


where M e = U e /c is the Mach number of the mean flow, which has been assumed small. 
The anti-symmetric component of the local flow past the leading edge lags the incident 
field (and the symmetric component) by a phase angle 7r/4, owing to the influence of 
compressibility. The anti-symmetric component again dominates the symmetric compo- 
nent, except for values of 9 near zero. However, the dependence on the incidence angle 9 
is fundamentally different than for the case of small k. The anti-symmetric component of 
the slip velocity has a peak value 2e' 7r ^ 4 (c/7ru>r n ) 2 sin kQ at the nose, a result that depends 
only on the compressibility of the fluid and not on the airfoil semi-chord b. 

In the next section, the symmetric and anti-symmetric components, C„(S) and C a (S), 
of the Receptivity Coefficient are calculated, and sample results are presented illustrating 
the dependence of the Receptivity Coefficient Ci(S) on the incidence angle 9 of the 
acoustic wave in the limits k <£. 1 and k 1. 


2.5 Numerical Results 

In §2.3, we utilized asymptotic methods to obtain the non-Blasius generalizations of the 
Lam-Rott asymptotic eigenfunctions, the first of which is the precursor of the Tollmien- 
Schlichting wave. The asymptotic analysis determines the form of these eigenfunctions, 
but not their coefficients C*. It appears that the Receptivity Coefficient C\ for a particular 
free-stream disturbance can be determined only by numerically solving the LUBLE over 
the full range of £ and examining the behavior for large £. 
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Figure 2.3: Plot of Q(£) = (Inf 2 ) l £ 2 ([4>" — F")/F") rt -o against (Inf 2 ) *. Comparison 
with the large f asymptotic form, (2.23) yields B\ « 2.08. 

In contrast to the flat-plate studies of previous authors, where the mean flow is known 
throughout, in the present study the mean flow, governed by (2.19), must be computed 
as well as the unsteady disturbance, governed by (2.21). The initial conditions at £ = 0, 
which correspond to the steady Hiemenz flow and its quasi-steady linear perturbation, 
take the form of ODEs in r /, which were solved using a fourth-order Runge-Kutta method. 
The solution was then obtained by marching downstream using a Keller Box scheme 
(Keller & Cebeci, 1970) for both the mean flow, given by a nonlinear PDE (2.19), and the 
linearized disturbance equations for symmetric and anti-symmetric components, (2.53). 
From the asymptotic form of the eigensolutions, it is clear that the growth or decay of the 
disturbance is very sensitive to the mean flow. For this reason, the departure of the mean 
flow away from the Blasius solution was calculated rather than the mean flow directly. 
The value of the unknown coefficient B x in (2.23) can then be extracted from the wall 
shear. From the asymptotic expansion for the mean flow close to the wall (2.25), 

(2 - 63) 

From figure 2.3, where the above quantity is plotted for real f, comparison of the 
numerical result with this asymptotic form gives B x « 2.08. This is somewhat higher 
than the value suggested by Van Dyke (1964b), who estimated B x « 1.6 - 1.9 based on 
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coarse numerical results and on numerical patching of the asymptotic series valid near 
the nose to the series valid far downstream. The computations undertaken in the present 
work extend to much higher values of f, allowing much more accurate extrapolation to fix 
B\. In addition, the magnitude of the next term in the asymptotic series can be estimated 
from figure 2.3. By considering the magnitude of the correction term at £ 2 = 30, the 
coefficient of the next term in expansion (2.63) is estimated to be approximately 10, 
which probably explains the under-estimation of Van Dyke. 

We turn now to the solution for the disturbance. Since the real part of T(£) is 
negative, the first Lam-Rott eigensolution becomes exponentially small far downstream, 
compared to the Stokes wave. Moreover, since the eigenvalues are inversely ordered, the 
first eigensolution is also exponentially small compared to all the other eigensolutions. 
Thus, it is very difficult to extract the coefficient of the eigensolution by direct numerical 
solution of the LUBLE. A way round this difficulty is to move the integration off the real 
line into the complex £-plane in such a way that the eigensolution will grow exponentially 
(Goldstein et al, 1983). This will occur if the streamwise variable £ is chosen such that 
— 57t/12 < arg £ < — 7t/12. The inverse ordering of the eigensolutions is also rectified 
by this process, the first eigenfunction becoming exponentially dominant for |£| 2> 1. 
‘Peeling off’ the exponentially growing part of the first eigensolution, that is solving for 
/ = x/je~ T i rather than xp, allows the numerical solution to be continued much further 
downstream while retaining accuracy (Heinrich &: Kerschen, 1989). 

There are essentially two ways of obtaining the Receptivity Coefficient C\, one based 
on the wall shear and the other based on the oscillating boundary-layer thickness. The 
latter method was used by Goldstein et al (1983) and Heinrich & Kerschen (1989) for the 
fiat-plate analysis. In the present work the two methods gave virtually identical results, 
small differences only arising due to the difference in large-£ extrapolation procedure in 
the two cases. Using the wall shear as the basis for comparison, we define 


C(S,Q = 01 , 

Jo 


(2.64) 


where = 0.4356(1 -I- i) (£ 2 /2) n+1 . Since the argument of £ was chosen so that the first 
eigensolution dominates other components of the solution, / Ciipie and hence 


C~ C 1 ! 


1 + b^ 1887 


c 0.774 


+ ^ + ^M23 +it 5>U... 


(2.65) 


as £ — ► oo. Thus in principle, C\ can be extrapolated by plotting the numerically obtained 
value of C(S, £) against £-° 774 . However, this method has serious limitations as far as 
obtaining highly accurate estimates for C\. Until f is extremely large, the sizes of the next 
three terms in the expansion are very close to that of the leading order correction, with 
the relative magnitudes depending also on S. In addition, the analysis of the base flow 
suggests that fq is quite large, so for larger values of S the exponential term cannot be 
expanded out until £ is very large. Due to these difficulties, the Receptivity Coefficient 
for non-zero S cannot be calculated to the same accuracy as for the flat plate, the 
inaccuracy increasing as S increases. In our calculations, we usually chose arg f = — 7 r /4 
which ensures that the rate of growth of the first Lam-Rott eigensolution is maximized. 
In addition, for this angle any error in the numerically determined value of B { enters 
only the phase of C\ and not the magnitude. For some values of S, the solution for 
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arg £ = — 7r/3 was also obtained and the extrapolated value of C\ was virtually identical, 
further verifying the numerical work as well as the functional form of the eigensolution. 
The solution was integrated up to |£| = 15, except for the larger values of 5 where it 
was continued up to |£| = 20 in an attempt to improve the extrapolation. Computation 
to larger values of £ becomes progressively more expensive due to the presence of the 
viscous wall layer, the relative thickness of which decreases downstream. 

In figure 2.4, the numerically obtained value of the Receptivity Coefficient for acoustic 
waves propagating parallel to the airfoil chord, C 31 is plotted for several values of 5. In 
figure 2.4a it is seen that the magnitude of the Receptivity Coefficient increases slightly 
for very small 5, but decreases rapidly as the nose radius increases further. For 5 = 0.3, 
the receptivity is reduced to approximately 15% of that for a flat plate. The decay 
of |C a | appears to be exponential in 5, though this has not been verified analytically. 
For 5 < 0.1, the results for \C 3 \ are estimated as accurate to 2%, while the accuracy 
for 5 = 0.3 is only about 10%, due to the difficulties discussed above. However, the 
large relative uncertainty for the latter case is of little practical importance, since the 
receptivity is so low by this stage. The small increase in the parallel-wave receptivity seen 
for very small values of 5 is in agreement with small-5 asymptotic theory which will be 
reported elsewhere (Hammerton & Kerschen 1996a). Figure 2.46 shows the change in the 
phase of the Receptivity Coefficient as 5 increases. The increase in arg(Cj) corresponds 
to a lag in phase of the instability wave. 

As we have already noted, computation must be carried out to large values of £ for 
accurate determination of the Receptivity Coefficient, since the eigensolution tends to 
its asymptotic limit slowly. However, the ratio of Receptivity Coefficients for different 
free-stream disturbances can be evaluated accurately at only moderately large values 
of £, since the same asymptotic behavior of the eigensolution arises whatever the form 
of the free-stream disturbance. In the present chapter, we consider only the effect of 
acoustic waves at oblique angles. In figure 2.5, the receptivity due to the symmetric 
and anti-symmetric components of the free-stream disturbance is compared for different 
nose radii. We define A = C a /C 3 , where C a is the Receptivity Coefficient related to the 
anti-symmetric component of (2.50). In figure 2.5a, |A| is plotted as a function of 5. 
We see that \C a \ is larger than |C,| by approximately a factor of six in the flat-plate 
limit 5 — + 0, where the anti-symmetric component of the slip velocity has a square-root 
singularity at the leading edge. The value of |A| decreases rapidly as a finite nose radius is 
introduced, relieving the singularity at the leading edge. Thus, the behavior of C a for very 
small 5 is quite different from that of C a , where a small rise in magnitude was seen. The 
sharp change in A for small 5 suggests the appearance of singular behavior in the small-5 
expansion, in contrast to the behavior seen in figure 2.4 for the parallel-wave case. Beyond 
5 = 0.05 the decline in |A| becomes more gradual, the value of |A| remaining above 2.5 
right out to 5 = 0.3. For all values of 5 examined, the anti-symmetric component of the 
slip velocity was found to be more effective than the symmetric component in generating 
a Tollmien-Schlichting wave. 

The phase difference between the C a and C s is plotted in figure 2.56. It is interesting 
to note that as 5 — > 0, arg A — >- 7t/4. For reasons given above, A can be calculated 
very accurately, and this result is true to an accuracy better than 0.1%. However, the 
fundamental reason for this result is not clear. Arg(A) initially drops rapidly as a finite 
nose radius is introduced, then remains roughly constant from 5 = 0.05 out to 5 = 0.3. 
The positive values of arg(A) correspond to a phase lag for the contribution from the 
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Figure 2 . 4 : Variation of the receptivity coefficient as a function of Strouhal number, 5 , 
for an acoustic wave in the free stream propagating parallel to the body axis: (a) the 
amplitude of the receptivity coefficient |C a (S)|; (b) the phase of the receptivity coefficient 
arg(C a )/7r. 
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Figure 2.5: Comparison of the relative receptivity due to symmetric and asymmetric 
components of an oblique acoustic wave, A = C a /C 3 : (a) the magnitude of the relative 
receptivity, |A|; (b) the phase, arg(A). 
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Figure 2.6: Variation of |Ci| with acoustic incidence angle 6 , for k 1 and non- 
dimensional airfoil chord length a = 10. The solid line is for S = 0.0, the dotted 
line S = 0.2. 


anti-symmetric component of slip velocity, relative to the contribution from the sym- 
metric component. This phase lag may be related to the fact that the anti-symmetric 
component of the slip velocity takes on its largest values near the nose, while the sym- 
metric component takes on its largest value farther downstream. The concentration of 
the anti-symmetric component near the leading edge is particularly pronounced in the 
limit S — ► 0, where arg(A) takes on its largest value, the slip velocity being singular at 
£ = 0 in this case. 

The variation of the total Receptivity Coefficient C\ with acoustic wave incidence 
angle is illustrated for the case k -C 1 in figure 2.6. A representative aerodynamic 
reduced frequency, a = 10, has been chosen, and results are plotted for two nose radii, 
5 = 0 and 0.2. The receptivity for the flat-plate case ( S — 0) is dominated by the 
contribution from the anti-symmetric component, causing the shape of the plot for |Ci| 
to be quite close to sin#, except in the vicinity of 6 = 0° and 180° where the level is 
determined by the symmetric component. The case a = 10, S = 0.2 corresponds to a 
typical airfoil design. The overall receptivity level for S = 0.2 is smaller than for the flat 
plate, due to decreases in both \C 3 \ and \C a \. Since the finite nose radius causes a larger 
decrease in C s than C a , the influence of the symmetric component of the free-stream 
disturbance is somewhat larger in this case. The peak receptivity for S = 0.2 is roughly 
one-fifth that for the flat-plate. 

The variation of the total Receptivity Coefficient C\ with incidence angle is illustrated 
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Figure 2.7: Variation of |Ci| with acoustic incidence angle 9, for k » 1 and M = 0.1, 
and the same values of S plotted in figure 2.6. 

for the case k 1 in figure 2.7. A representative Mach number M e = 0.1, has been 
chosen. The dependence of the receptivity level on acoustic wave incidence angle is 
very different from the case k < 1 illustrated in figure 2.6. The flat-plate result is 
again dominated by the contribution from the anti-symmetric component of the free- 
stream disturbance, but the plot-shape resembles sin ^9, except near 9 = 0° where the 
symmetric component determines the level. As before, for S = 0.2 the overall receptivity 
level is decreased, but the relative contribution from the symmetric component of the 
disturbance is somewhat larger. 

2.6 Conclusion 

Boundary layer transition is influenced both by the stability properties of the boundary 
layer, and by the characteristics of the free-stream disturbances and the receptivity of 
the boundary layer to these disturbances. The present work has examined the influence 
on leading-edge receptivity of the nose radius of an uncambered airfoil symmetric mean 
flow, in the low Mach number limit. 

Leading-edge receptivity involves a gradual evolution of the boundary-layer distur- 
bances with downstream distance, eventually leading to a growing Tollmien-Schlichting 
wave. The Receptivity Coefficient is essentially the amplitude of the asymptotic eigen- 
function which is the precursor of the Tollmien-Schlichting wave. It is the Receptivity 
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Coefficient which contains all the relevant information about the free-stream disturbance. 
The decay rate of the asymptotic eigenfunction, like that of the Tollmien-Schlichting 
wave, is influenced only by the mean flow. One of the most important properties is the 
pressure gradient parameter, j3. For symmetric flow past a parabola, 

e=?TS’ ( 2 . 66 ) 

where S is the Strouhal number, S = u>r n /U e . Thus, the pressure gradient is everywhere 
favorable, decreasing monotonically from its maximum value at the stagnation point and 
approaching zero far downstream. Well downstream of the nose, the favorable pressure 
gradient has a stabilizing influence on the boundary layer, leading to increased damping 
of the Tollmien-Schlichting waves and of the precursor disturbances. In addition, the 
stronger pressure gradient close to the nose has an effect on the Receptivity Coefficient. 
As the nose radius of the body is increased (i.e. S increased), the strongly favorable pres- 
sure gradient near the nose extends over a larger number of disturbance wavelengths, and 
this additional stabilizing influence is likely to result in a decrease in the Receptivity Coef- 
ficient. Numerical calculations show that the symmetric component C 3 of the Receptivity 
Coefficient decreases rapidly with increasing S, dropping to approximately 12% of the 
flat-plate value when S = 0.3. The anti-symmetric component C a decreases even more 
rapidly, dropping to 5% of flat-plate value when S = 0.3. The total Receptivity Coeffi- 
cient Ci is much larger for oblique acoustic waves than for parallel acoustic waves, since 
the coefficient K a (6) that multiplies the anti-symmetric component C a is large compared 
to the coefficient k 3 (9) that multiplies the anti-symmetric component C 3 . 

While the detailed asymptotic analysis depends on the exact geometry of the leading 
edge, the rapid decrease in receptivity with increasing nose radius should also be valid 
for other streamlined bodies. However, symmetric mean flows past leading edges less 
streamlined often involve regions of adverse as well as favorable pressure gradient close 
to the nose. For such bodies, it is possible that an increase in leading edge thickness 
could lead to an increase in the Receptivity Coefficient. For an airfoil at an angle of 
attack, both favorable and adverse pressure gradients exist in the neighborhood of the 
nose. Analysis of the receptivity in this case will be presented in chapter 3. 

Comparing the asymptotic analysis presented here with numerical and experimental 
results is difficult. Numerical work by Murdock (1981) for a parabola in a flow with 
a parallel acoustic wave also showed a decrease in receptivity as the nose radius was 
increased. The computations of Lin et al. (1992) were for parallel acoustic waves incident 
on half-ellipse leading edges connected to a flat plate and for super-ellipse leading edges 
(which avoid the discontinuity in curvature), geometries chosen to match the experiments 
of Saric et al. (1994). For both these geometries there are regions of adverse pressure 
gradient near the leading edge, as well as the possibility of additional localized receptivity 
mechanisms (Goldstein, 1985); hence no direct comparisons can be made with the present 
results. In addition, the receptivity level ‘seen’ in experiments and full Navier-Stokes 
calculations is a combination of the receptivity process described in the current chapter, 
and the stability characteristics farther downstream. For the flat-plate, Goldstein (1983) 
provided the asymptotic analysis of the triple-deck region and thus (in theory) could 
calculate the disturbance amplitude at the lower branch. Reproducing such an analysis 
for a general airfoil surface would provide little additional physical insight concerning 
the receptivity process. By restricting attention to the region in which forcing by the 
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free-stream disturbance occurs, the process of receptivity is isolated from the instability 
phenomena that occur farther downstream. For a global picture of the transition process, 
the current analysis provides an upstream boundary condition for linear stability analyses. 
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Chapter 3 


Low frequency theory for a leading 
edge with no aerodynamic loading 

3.1 Introduction 

The receptivity process through which free-stream disturbances generate instability waves 
in boundary layers generally comes about through non-parallel mean flow effects, which 
may arise either in the leading-edge region, or in a localized region farther downstream 
in the boundary layer (Goldstein & Hultgren 1989; Kerschen 1990). 

In Goldstein (1983), an asymptotic analysis was developed for leading-edge recep- 
tivity on a semi-infinite zero-thickness plate. Leading-edge receptivity coefficients for 
various free-stream disturbances were calculated by Goldstein, Sockol & Sanz (1983) and 
Heinrich & Kerschen (1989). However, aerodynamic bodies designed for subsonic flow 
generally have finite thickness distributions with a parabolic leading edge. In chapter 
2, the influence of the thickness of a body was examined. The situation considered was 
a thin, symmetric airfoil of chord 2b at zero angle-of-attack in a uniform flow of speed 
U . A plane acoustic wave of frequency u, propagating at an angle 6 with respect to 
the airfoil chord, was assumed to be incident on the airfoil as illustrated in figure 3.1. 
A high Reynolds number asymptotic analysis (e 6 = uu/U 2 \ e < 1) was formulated for 
an incompressible, two-dimensional flow. The nose radius of the leading edge, r„, enters 
through a Strouhal number, 


which is the ratio of the nose radius and the hydrodynamic length scale, U/u. 

When S = 0(1), two streamwise regions enter the analysis, one region where the dis- 
tance downstream is 0{U /u) and the disturbance is governed by the linearized unsteady 
boundary layer equation (LUBLE), and a second region at distances 0(e~ 2 U/u) where 
the disturbance is governed by the triple-deck structure, corresponding to the asymptotic 
form of the Orr-Sommerfeld equation (OSE) in the vicinity of the lower branch. In the 
LUBLE region, the inviscid pressure field and slip velocity induced by the free-stream 
disturbance drives the unsteady motion in the boundary layer. Far downstream in the 
LUBLE region, the solution consists of a generalized Stokes wave, and a set of asymp- 
totic eigensolutions. The wavelengths of these eigenfunctions shorten progressively with 
distance downstream. Eventually, the self-induced pressure field associated with the 
displacement thickness of each asymptotic eigenfunction becomes significant, and the 
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Figure 3.1: An illustration of the physical situation of interest: a thin, symmetric airfoil 
of chord 26 is at zero angle-of-attack in a uniform flow of speed U, with a plane acoustic 
wave incident at an angle 6 with respect to the airfoil chord. 
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triple-deck structure replaces the LUBLE as the correct asymptotic approximation to 
the Navier-Stokes equation. It can be shown that the first asymptotic eigenfunction 
of the LUBLE matches on to the Tollmien-Schlichting wave solution of this triple-deck 
region. Thus, the form of the free-stream disturbance and the geometry close to the nose 
influence the amplitude of the Tollmien-Schlichting wave only through the coefficient C\ 
of the first asymptotic eigenfunction. Therefore, we call C\ the ‘receptivity coefficient’. 
The asymptotic analysis far downstream in the LUBLE region determined the form of 
the eigenfunctions, but not their coefficients C*. These coefficients can be found only 
through a full solution of the LUBLE, which must be determined by numerical methods. 

Theory for 5 = 0(1) was developed in chapter 2, where associated numerical results 
were also presented. These results showed that the receptivity associated with a sym- 
metric free-stream disturbance increases slightly with 5, for small 5 (figure 2.4a), but 
the receptivity associated with a free-stream disturbance that is anti-symmetric about 
the leading edge decreases sharply with 5 (figure 2.5a). In the present chapter, a small-5 
asymptotic theory is developed. In this limit, the streamwise development of the bound- 
ary layer can be divided into three regimes. In the nose region, where the distance 
downstream is 0(r n ), the unsteady motion in the boundary layer is quasi-steady at lead- 
ing order. In the LUBLE region, where the downstream distance is 0(U/u ), the unsteady 
terms enter at leading order in the disturbance equations, but the mean pressure gradi- 
ent is small so that the mean boundary layer flow approaches the flat-plate solution. In 
the OSE region further still downstream, a triple-deck structure arises as before. In the 
small-5 asymptotic analysis presented in this chapter, we focus on the nose and LUBLE 
regions. The subsequent asymptotic matching to the OSE region follows that discussed 
in chapter 2. Our analysis involves two small parameters, 5 and e. Formally, we take 
6 — » 0 and then consider the small-5 limit, obtaining correction terms of 0(5 1 / 2 , 5). We 
require the boundary-layer approximation to be valid in the nose region xj = 0(r n ), 
which is true if the Reynolds number based on nose radius is large, 


Ur n 

V 



(3,2) 


This condition is also equivalent to the requirement that the perturbation in the recep- 
tivity coefficient due to finite Reynolds number is smaller than that due to geometric 
effects. Thus the asymptotic theory is more relevant to the limit ui — ♦ 0 rather than 

r n -* 0. 

In §3.2, the equations governing the mean flow and the time dependent perturbation 
in the boundary layer axe obtained. The small-5 asymptotic structure of the unsteady 
flow in the boundary layer, produced by the symmetric and anti-symmetric components of 
the free-stream disturbance, is analyzed in §§3.3 and 3.4 respectively. In §3.5, numerical 
solutions of the asymptotic equations are presented, and the resulting asymptotic expres- 
sions for the receptivity coefficient are compared with the numerical results of chapter 2. 
Finally in §3.6, the results are summarized and brief comparisons with experiments are 
made. 


3.2 Formulation 

Full details of the derivation of the equations governing the boundary layer flow are given 
in chapter 2, together with interpretation of the evolution of the unsteady disturbances 
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that eventually lead to instability. Here we provide only the key equations necessary to 
illustrate the small Strouhal number limit which is the subject of the present chapter. 

The streamfunction ^ in the boundary layer, non-dimensionalized by the quantity 
U is 2 is expressed in the form 

9 = £<P + r/fe~ il (3.3) 

where the unsteady component ip is assumed small compared to the mean flow component 
£<p. We use parabolic coordinates (£, 77) defined by 

x d + Wd = " ([^ + i(S l/2 + e 3 77)] 2 + s'j , (3.4) 

where x<f, yd are (dimensional) Cartesian coordinates centered on the airfoil leading edge. 
The parabolic coordinate £ along the airfoil surface has been non-dimensionalized by the 
square root of the disturbance length U/u. For the parabolic coordinate 77 in the normal 
direction, the origin has been shifted to the airfoil surface, and the additional factor e 3 in 
(3.4) corresponds to non-dimensionalized by the quantity {v/U)i. Thus 77 corresponds 
to the conventional normal coordinate for the mean boundary layer. 

3.2.1 Steady flow 

The steady boundary-layer flow is most naturally expressed in terms of a streamwise 


variable scaled on the nose radius, 

Z = Z/S l/2 . (3.5) 

The function <£(£, 77) describing the steady boundary-layer flow then satisfies the differ- 
ential equation 

<P + £(07^$ — — j — 1) = Oi (3-6) 

with boundary conditions 

(p = (pr, — 0 at 77 = 0 and <p v — ► 1 exponentially as rj — > oo. (3.7) 
For large f, <p takes the form (Van Dyke, 1964) 

<Kl V) ~ F{ 77) + G^H- + G 2 (77)i + O , 72 ~ 3.774, (3.8) 


where the order of the next higher order term is discussed below. In this expansion, F(tj) 
is the Blasius function, and the functions G* (77) satisfy homogeneous boundary conditions 

Gj(0) = G((0) = 0, G[ — ► 0 exponentially as 77 — * 00, (3.9) 

with governing equations 

C? 2 (Gi) = 0, G 2 (G 2 ) = F ' 2 - 1 + 2 (F"G, - F'G\), (3.10) 
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where we define the set of operators by 

QiU) = /'" + Ff" + 7 F'f' + (1 - T )F7. (3.11) 

In the first perturbation term of the large-£ expansion we have Gy — AyGy, where 

G\ = t)F'(t]) - F(r]), (3.12) 

which arises as an eigenfunction of Q 1 for 7 = 2. The coefficient Ay remains undetermined 
at this order. However, the equation for G 2 (t]) then takes the form 

g 2 (G 2 ) = F' 2 - 1 + 2 A l FF". ( 3 . 13 ) 

Since the Blasius function F is the solution to the adjoint of Q 2 , the value of Ay is fixed 
by the solvability condition for G 2 (tj) which leads to 

f°° F(F' 2 - 1 + 2A l FF")dr ] = 0, (3.14) 

J 0 

and gives Ai = 0.60115. 

The equation for G^tj) involves the same operator Q 2 and hence contains the same 
eigenfunction. Thus, G 2 = ByGy 4- G 2 , where G 2 (rj) satisfies (3.13), with boundary 
conditions (3.9) and G^O) = 0. Since £~ 2 G 1 is an eigensolution of the boundary-layer 
perturbation equation, the value of the coefficient By appearing in the 0(£~ 2 ) term can 
not be determined by the large-£ analysis and therefore depends on conditions close to the 
nose of the body. Numerical integration of (3.6) from £ = 0 gives By « 2.08 (chapter 2). 
For all 7 > 0, Q 7 has solutions which decay algebraically as 77 — ► 00, but only for certain 
7 do solutions exist which decay exponentially at infinity. Libby h Fox (1963) give the 
first 10 eigenvalues. The first such eigenvalue, 71 = 2 has already been discussed. The 
next eigenvalue, 72 ~ 3.774, leads to the next higher order term in the large-f expansion 
(3.8). 


3.2.2 Unsteady flow 

The tune-dependent contribution to the boundary-layer flow satisfies the linearized un- 
steady boundary layer equation (LUBLE), 


^) = (5 + e 2 ) 1/2 


i(s+e)~ 


s + C 


F(ip) = Vv™ + [0 + + 

with boundary conditions 


\ du 3 ' 

;] u ” 


i(S + (‘) - - 


■s-e 
s + Z 


2 <t>n 


fpr, 


(3.15) 


ip = 'tp n = 0 at 7] = 0 and ip n —> (S + ^ 2 ) 1/2 u 3 as rj -* 00 . (3.16) 

Here, u s (£) is the slip velocity induced on the outer edge of the boundary layer by the 
free-stream disturbance. Note that this equation is written in terms of £, the streamwise 
coordinate scaled on the disturbance length scale U/u, and hence (j> must be expressed 
in terms of f. 
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The development of the solution ip is described in detail in chapter 2. Far downstream 
(£ ^ l)i ip consists of a particular solution, xp p , driven by the local value of the unsteady 
pressure gradient, together with an infinite set of asymptotic eigensolutions, 

S) = ^( 17 , ft 5) 4- £ Ci(SMfo,ft 5). (3.17) 


The eigensolutions, ipi, depend on the geometry of the body far downstream, but are 
independent of the local free-stream disturbance. In this chapter, where we are concerned 
with the generation of instability waves in the boundary layer, we consider eigensolutions 
which are generalized forms of the Lam &; Rott (1960) eigenfunctions. Another set of 
eigensolutions could also be calculated, corresponding to generalizations of the functions 
derived by Brown & Stewartson (1973) for a flat plate. The role of these two sets of 
eigensolutions is discussed in chapter 2. The coefficients Cj multiplying the asymptotic 
eigenfunctions are determined entirely by conditions close to the leading edge (£ = 0(1)). 
One of these eigensolutions, which we label ipi, matches on to the Tollmien-Schlichting 
wave in the Orr-Sommerfeld region farther downstream, where £ = 0(e -1 ). Thus, it 
is only through the coefficient C\ that the unsteady disturbances in the free stream 
influence the amplitude of the Tollmien-Schlichting wave. Our primary interest is in the 
relationship between the free-stream disturbances and the amplitude of the Tollmien- 
Schlichting wave, as a function of Strouhal number S. Thus, we focus on the asymptotic 
eigensolutions of the LUBLE, which develop a two-layer structure for £ 1. From 

chapter 2, equations (2.31), (2.48) and (2.49), at the outer edge of the boundary layer 
(77 -+ 00), the first asymptotic eigenfunction takes the form 


rf 0) = -0.69213, rf 1 ' 

T}°\0 = -A? 1 ji - 


(i+o(r°™)), 

= -1.9878 i, 

+ (2Aj + f - Bi) — 


\ 



(3.18) 


where A = /{(?\ 2 U’ 0 ). Here pi « 1.0187, the first root of Ai'(— p) = 0, and U' Q = 

F"(0) ~ 0.4696. The numerical constants A\ and B\ arise from the asymptotic form of 
the mean boundary- layer flow far downstream, (3.8). 

For a low Mach number flow, outside the boundary layer the unsteady flow in the 
vicinity of the leading edge is incompressible and irrotational. Potential flow theory then 
shows that this local flow consists of symmetric and antisymmetric components of the 
form 

U *(0 = K 3 + £ 2 ) 1/2 + Ka (5 + £ 2 ) 1 / 2 - ( 3 - 19 ) 

Here ^/(-^-l-^ 2 ) 1 ^ 2 and 1/(5 + £ 2 ) 1/f2 correspond to purely symmetric and anti-symmetric 
flow about the leading edge, respectively. The coefficients k s (9) and K a (&) multiplying 
the symmetric and antisymmetric components are independent of the nose geometry, 
but depend on the free-stream disturbance, being determined by global features of the 
unsteady flow. When the unsteady disturbance velocity is parallel to the mean flow, (i.e. 
0 = 0), the forcing is entirely symmetric and k 3 = 0. For oblique (0 ^ 0) disturbances, 
the unsteady flow will have both symmetric and anti-symmetric components. Calculation 
of the coefficients k, and K a for different free-stream disturbances is discussed in §2.4. 
Writing ip = n,ip s + K a tp a , and substituting into (3.15), it follows that the receptivity 
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coefficient can be decomposed into contributions from the symmetric and anti-symmetric 
components of the unsteady outer flow, 


Oi(5) = k s C 3 (S) + K a C a (S), 

(3.20) 

where C, and C a are obtained from the solutions of 



(3.21a) 


(3.216) 


respectively. The numerical methods used to obtain values of C a (S ) and C 3 (S) for fixed 
S are described in chapter 2. The remainder of this chapter is concerned with the 
asymptotic expansions of C a (S ) and C 3 (S ) as S — ► 0. 

In the nose region, f = 0(1), the mean pressure gradient decreases with distance 
downstream, and the mean flow approaches the flat-plate limit as £ — ► oo. In the 
small-5 limit, the nose-region is quasi-steady at leading order; the linearized boundary 
layer equations only become fully unsteady when f = 0(1), i.e. £ = 0(S~ 2). Thus 
in the receptivity region f = 0(1), the mean flow is close to the Blasius solution, and 
the receptivity coefficient is calculated as a perturbation away from the flat-plate value. 
Receptivity to the symmetric and anti-symmetric components of the unsteady free-stream 
disturbance is analyzed in §§3.3 and 3.4, respectively. In each case the solutions in the 
nose and receptivity regions are asymptotically matched, and then the large-f behavior in 
the receptivity region is compared to (3.18) in order to extract the receptivity coefficient, 
C 3 or C a . 


3.3 Receptivity to symmetric forcing 


For symmetric forcing, (3.21a) suggests that for small S the unsteady perturbation in 
the nose region, £ = 0(1), is quasi-steady at leading order and takes the form 


= sH (e Q & n) + 5 6,(1 t}) + o(s 2 )) , 

where the functions 6 0 and 9 X satisfy 


M 3 (6 0 ) 
M 3 (6,) 


1 + £ 2 ’ 

K 1 +0(1- (W, 


and 


(3.22) 


(3.23) 


M s (p) = Prjw + (0 + £,<t>i)Prm + Pn + + ZihmPi ~ ( 3 ‘ 24 ) 

The factor is extracted from tp 3 in (3.22) to simplify the boundary conditions, which 
become 


0 

II 

0 9e ‘ 

d6 0 

drj 

36, 

dp 


on 77 = 0 
exponentially as 77 - 


(3.25) 


00. 
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Here we are interested in the solution in the large-£ limit, in order to provide upstream 
boundary conditions for the receptivity region f = 0(1). In this limit, the mean flow is 
given by (3.8) and the unsteady perturbation takes the form, 


Bo ~ H 0 (v)+ l ^H 3 ( V ) + lH i (T 1 ) + 0(C y2 ), 1 
0 X ~ + h^f 2 )//^) + H 4 (t)) + 0(^ 2-72) ), J 


(3.26) 


where we have adopted this rather unusual labeling notation in order to retain consistency 
with the notation used in §3.4. The set of functions (77) and Hi(j]) satisfy 

Go{H 0 ) = 0, G- 2 (Ho) = 

G 2 (H 3 ) = ttg(J5To), G 0 (H 3 ) = 

g 2 (H 4 ) = -2-n^(H 0 ) Go(H 4 ) = 

+2 F'H' 3 - 2F"H 3 , 

subject to boundary conditions 

H t = H'= 0, Hi = H'i = 0, 
i/'-l, H’^-* 0 , //'- 0, 

where the decay is exponential as 77 — » 00. Here the operator £? 7 is defined in (3.11) 
and the operators 7Z^ are defined in Appendix C. As was noted in §3.2, Q 2 has an 
eigensolution (i.e. a solution of G 2 (H) = 0 with homogeneous boundary conditions). 
Thus H 3 and H 4 can not be determined entirely by the large-£ asymptotic analysis, but 
require knowledge of the full solution in the nose region. However, for the symmetric 
case considered in this section, the leading-order term in the small-5 expansion can be 
determined exactly in terms of the steady flow (Lighthill 1954). In the notation of the 
present chapter this gives the quasi-steady solution as 0 O = |(^0, + 0) and hence, 

Hq — zivF + F), H 3 = + (?i), H 4 = 5(77^2 + G 2 ), (3.29) 




>(i - K), 

-iH' 3 +-Rgl(H 0 ), 

+2 F'H's - r 2F"H s , 


(3.27) 


on 7] — 0, 
as 77 —► 00, 


(3.28) 


where Gi(tj), G 2 {rj) are defined in (3.12, 3.13). Solutions for Hi are obtained using a 
fourth-order Runge-Kutta method, shooting from 77 = 0 and using the requirement of 
exponential decay as 77 — » 00. 

It is clear from (3.26) that the small-5 expansion (3.22) breaks down when 5£ 2 = 
0(1), which is to be expected since when £ = 0(1) we have reached the receptivity 
regime where unsteadiness enters at leading order. Within the receptivity region, we are 
interested only in the coefficient of the first eigensolution. We therefore anticipate the 
large-£ form (3.18) by setting ip = £u>(£, r7)e _Af3/3 , where A = e~4 ™ /{p\ ,2 Uq) as before. 
From (3.21a), w is given by 


M s {w) = 


Ks + t 2 )- 


K(w) — Wjjqri -f- [0 + + 


25 

5 + £ 2 


,A4 3 /3 


i(S + ?)-&*+ [K 3 - 


25 


+ (1 - A ^)(prjr,W + - 0„UV;). 


5 + C 2 


07) 


W r 


(3.30) 
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Using the large-£ form of <j> (3.8), we see that the coefficients which appear in Af 3 ( w) 
involve terms in (£ 2 /5)~ l ln(£ 2 /5) and {£ 2 /S)~ l , and hence we can expand M s as an 
asymptotic series in 5 , 

K(w) ~ A ri°!(w) + 5 In S^Siw) + SMgi(w) + 0(S” /2 ), (3.31 ) 

where the partial differential operators Nq^{w) are defined in Appendix C. Now writing 

to(f, 77 ; 5) ~ wo + 5 In Sw 3 + Sw 4 + 0{S’ n/2 ), (3.32) 

and expanding the right side of (3.30), we obtain a set of equations governing the evolution 

Of Wi(£,T)), 

Klim) = 

KiM = -A® (i»o), ' 

<?(«.) = (i-2 

Matching back to (3.26) we see that as f — ► 0, 

w 0 ~ Ho + £ 2 Hq, 

W3 ~ -75H3 - H 3 , 

w 4 ~ l -^H 3 + ±H 4 + \n(e)H 3 + H 4 . 

Thus we write 

Mo = e« J /3 + 

w 3 = — ^ H3 + < 73 , 

e 1 I 

w 4 = — ln(£ 2 )u > 3 + — e A? 3 / 3 /f 4 + q 4 , 

f 2 

and it can be shown after some algebraic manipulation that 

= ie* 3 ' 3 (l - #», 

<’(*>) = 

<?ta) = ie A ^ : >(i-J/;-H;)-i«Vo-X < 0 %) + 2F" g 3-2F' 9 ;, J 

where < 7 ' denotes the partial derivative of < 7 ,(£, 77 ) with respect to 77 . From (3.34), the 
functions qi satisfy initial conditions 

90 ( 0 , 77 ) = H 0 , 93 ( 0 , 77 ) = -H 3 (t)), 94 ( 0 , 77 ) = 7 / 4 ( 77 ), (3.37) 

and homogeneous boundary conditions, < 7 ,(£, 0 ) = ^(f, 0 ) = 0 , < 7 ' — » 0 exponentially 
as 77 —> 00 . The differential operators A/^, related to the operators A/^, are defined 
in Appendix C. Thus the evolution of the 0(5° ), 0(5 In 5) and 0(5) components of 
the unsteady flow has been cast in a form suitable for accurate numerical computation. 
Numerical solutions of (3.36) are discussed in §3.5. 


(3.33) 


(3.34) 


(3.35) 


(3.36) 
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3.4 Receptivity to antisymmetric forcing 

At first sight, the analysis of the receptivity in the anti-symmetric case appears to follow 
very similarly to the symmetric case. However, a subtle difference is seen in the structure 
of the solution. Motivated by the anti-symmetric forcing (3.216), we write the small-5 
expansion of the unsteady perturbation in the nose region, £ = 0(1), as 

*Pa = 4>o(£,v) + S {£,v) + 0(S 2 ). (3.38) 

The perturbation is again quasi-steady at leading order in the nose region, with 


•M„(« = C-±, 

M a ty l) = i(l + C 2 )(l — tyo)r)), 


(3.39) 


where 

M a (p) = M s {p) + (frrjP,, - 4> m p. (3.40) 

However, in contrast to the symmetric case, the large-£ limit for the unsteady flow 
now takes the form 


* ~ Mv) + Mn ) + ~Mr>) + ^Mv) + ~Mn) + 

4>i ~ £*Jo(v) +On(£ i )Ji( 1 l) +£J2 (v) + lnff 2 ) J 3 (t/) + J A ty) + 0(^ (3-72) ) 
where the functions J t ( rj) and Jity) satisfy 


(3.41) 


Gi(J a ) = 1) 

G2V1) = o, 

g 2 (J 2 ) = 2F'J[ — 2F"J l , 

Wz) = fcgOfr), 

QzW = -2-K§>(J 0 ) 

+2 F'4 - 2F" J 3 , 

with boundary conditions 


G-xtyo) = i(l - JS), 

Go(J i) = -iJ[, 

Gq(J 2 ) = -iJ' + 2F'J[ - 2F"J U 

Gxtyz) = + 7Z^l(Jo), 

Gi(J 4 ) = i(i - J' - J') - n^(j 0 ) 
+2 F'j' 3 - 2F"J 3 , 


(3.42) 


Ji = J[ = 0 Ji = J- = 0 on T] = 0 
J'q —* 1 , J/>! — > 0 J[ — > 0 as 77 — ► 00 

Compared to the expansion for the symmetric case (3.26), the additional terms in ln(£ 2 )/£ 
and l/£ are included in the large-£ expansion of xpo because the operator G 2 possesses 
an eigensolution. Thus J\ty) = PxtyF' — F), where the numerical constant Pi can be 
determined from the solvability condition for J 2 . As noted in §3.2, the Blasius function 
F is the solution to the adjoint of G 2 and we see that 

roc roc 

0 = / F(F'J[ - F"Ji) drj = 2Pi [ F 2 F" d 77 
J 0 Jo 

= Pi. (3.44) 

Hence J y = 0 and it follows that J x = 0, and J 2 (rj) = P 2 {t}F' - F), where P 2 is a 
numerical constant. Solutions to the set of equations (3.42) are obtained by shooting, as 
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described in §3.3, but the value of the constant P 2 in the expression for J 2 can only be 
determined by comparison with numerical solutions of ipo obtained by integrating (3.39) 
forward from £ = 0. Comparing the computed wall shear ^o(£-0), with the asymptotic 
form (3.41) gives 

P 2 ~ 6.07, (3.45) 

but accurate extrapolation is difficult due to the presence of higher-order terms in the 
expansion (3.41). The numerical results suggest that the coefficients of these terms are 
in fact relatively large. This is discussed further in §3.5. 

In the receptivity region £ = 0(1) where unsteadiness enters at leading order, we 
again anticipate the large-£ form of the eigensolution by setting ip a = ue -A ^ 3 / 3 . The 
function v(£, 77; S ) then satisfies 


K( v) = 


[ 1(s + f2 >-s4] 

e A P/ 3 
6 ’ > 

= a r.{v) + - 4> 

'rm V i . 


(3.46) 


where Af s (v) was defined in §3.3. As for the symmetric case, we expand the operator as 
an asymptotic series in S, 


K(v) + S\nSAfjjg(v) + SM$$(v) + 0(S^ /2 ), (3.47) 


but in this case, guided by the asymptotic form (3.41) in the nose region, we expand 
v{£,t}]S) as 

v ~ v 0 + S 1/2 v 2 + 5 In Sv 3 + Sv 4 + 0(S (72 " 1)/2 ). (3.48) 

Matching back to the nose region then suggests that we set 


Vo = e A?3/3 J 0 + £ 2 po. 

V2 = ^e^ 3/3 J 2 + fr?2, 

^3 = ^ 3 ^3 +P3, 

Vi = - ln(C 2 )u3 + -^e A * 3/3 Ji + p 4 . 

This finally leads to the set of equations 

-A/l°u(Po) = ie A?3/3 (l - Jo), 

No ;?(P2) = -ie^ 3/3 J^, 

^i°o(P3) = ie^J'-ATi^Po), 

N§(p 4 ) = ie A * 3 / 3 (l- -N { $(p 0 ) + 2F"p 3 -2F%, , 


(3.49) 


(3.50) 


with initial conditions 


Po(0,t}) = J 0 , p 2 (0, 77) = J 2 , P3(0, 77) = -J 3 (t7), p 4 (0,T]) = J 4 (tj), (3.51) 

and the homogeneous boundary conditions, p ; (£,0) = p'(£, 0) = 0, p\ — + 0 exponentially 
as 77 —> 00. As in §3.3, we use p' to denote the partial derivative of p,(£, 77) with respect 
to 77. Numerical solutions to this set of equations are discussed in the next section. 
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3.5 Numerical Results 


Having obtained expressions for the evolution of the unsteady flow in the form of an 
asymptotic expansion for small S, we now carry out the numerical integration required 
to obtain the corresponding asymptotic expansion for the receptivity coefficient Ci(S). 
As explained in chapter 2, in order to reliably extract the coefficient C\(S) of the first 
asymptotic eigensolution from the numerical results, it is necessary to carry out the inte- 
gration in the complex-^ plane, in the region -57t/12 < arg(£) < — 7r/12 where the first 
eigensolution of the LUBLE exhibits exponential growth and is the dominant component 
of the solution at large The sets of equations (3.36) and (3.50) were integrated forward 
in f, from £ = 0 to |£| = 10 using a Keller Box scheme (Keller & Cebeci (1970)), with 
arg(f) chosen as either — n/3 or — n/4 to ensure the dominance of the first eigensolution. 

For large £, the first asymptotic eigensolution is given by (3.18). Expanding the 
expression for the first eigensolution (3.18) for small S, and writing the receptivity coef- 
ficient, C a (S) or C 3 (S), as an expansion for small S , 

C(S) ~ C(0)(1 + c 2 S 1/2 + c 3 5 In 5 + c 4 S + ...), (3.52) 

we see that as f — * oo in the region of the complex plane of interest, 

0 ~ 0o (1 + C 2 S 1/2 + [C 3 — Aj4i£]*S In S + 

[c 4 + n ln(£ 2 /2) - \Z(2Ay + I-B 1 -A l ln(£ 2 ))]S), (3.53) 

where 0 O is the first eigensolution for the flat-plate case (Goldstein 1983). The coefficients 
C(0), c 2 , C 3 and c 4 can then be obtained by comparison with the numerical solution either 
at the outer edge of the boundary layer or at the wall. 

For the symmetric case clearly c£ s) = 0, while (7,(0), c[ 3 \ c 4 j) are extracted from the 
limiting forms 

(0) 

w o ~ C.( 0)(e/2) Tl [l + or 1 )], 

w 3 /w 0 ~ cf} 3 ) - A A^, > (3.54) 

Wa/wo ~ cJ' ) + Ti (1, ln(C 2 /2)-Ae(2i4 1 + |-B 1 -i4 1 ln(C 2 )), . 

as 77, £ — > 00. The evolution of w 0 is the flat-plate case described by Goldstein and gives 

C,(0) = -0.441 + 0.841 i (3.55) 

as the receptivity coefficient for a flat plate, in agreement with earlier results (Goldstein 
et al. 1983; Heinrich & Kerschen 1989). Using arg(f) = -7t/4, the corrections taking 
account of small nose radius effects are calculated to be 

c ( 3 s) « 0, c5 s) « 3. 12 - 2.37 i. (3.56) 

Integration using arg(£) = — 7r/3 gives a difference of less than 2% in the estimate for 
c{ 3) . Hence 

R e(C s (S)) ~ —0.441 -(- 0.62 S + 0(S°") 
lm{C 3 (S)) ~ 0.841 + 3.67 S + 0(S a '), a { = \ l2 ^ 1.887. 

In figure 3.2 asymptotic results (3.56) are compared with the results of the full integration 


(3.57) 
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Fig 2a — Real, symmetric Fig 2b — Imag, symmetric 



Fig 2c — Real, anti-symmetric Fig 2d — Imag, anti-symmetric 



Figure 3.2: Comparison of numerical and asymptotic results for the receptivity due 
to symmetric and anti- symmetric disturbances. Numerical results for the perturbation 
f(S) = (C(S) - C(0))/C(0) (see text) are denoted by o, asymptotic theory is marked by 
the solid line. Results for R e(/ s ) and Im(/ 3 ) are illustrated in (a) and (£>), respectively, 
and Re(/ a ) and Im(/ Q ) are illustrated in (c) and ( d ). 
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for small-5 presented in chapter 2. To make the comparisons as clear as possible, we plot 
the perturbation from the flat plate result f(S) = ( C(5 ) - C(0))/C(0) as a function of 
5 and compare this with the asymptotic form (3.52). In figure 3.2a and b, comparisons 
are made for the symmetric case, making it clear that the leading-order behavior for the 
symmetric receptivity coefficient is indeed described by (3.57). The small discrepancy 
between numerical results and asymptotic theory is likely to be due either to numerical 
uncertainty in extrapolating C(S) from large-£ calculations (as explained in chapter 2), 
or to the higher-order terms in (3.57). 

For the anti-symmetric case, 

v 0 ~ 

v 2 /vq ~ 
v 3 /v 0 ~ 
v 4 /v 0 ~ 

as 77, £ — i ► 00. Comparing these large-£ asymptotic forms with numerical results gives 

C a (0) « -5.33+ 1.66 i, (3.59) 

as the flat-plate solution, with small-5 corrections 

4 a) « -1.728 - 1.728 i, cj 1 * » 0, cj° « 3.12 - 3.57 i. (3.60) 

Again there is no significant difference in these results if arg(£) is varied. Hence 
Re(C a (5)) ~ -5.33 +12.08 5 1 / 2 -10.7 5 +0(5 Q2 ) } 

Im(G,(S)) ~ 1.66 + 6.34 S 1 ' 2 +24.2 S +0(S n ), a 2 = 1(72 - 1) « 1.387 ) 

(3.61) 

In figure 3.2c and d the asymptotic results (3.60) are compared with the results of the 
full integration for small-5 presented in chapter 2. It is clear that the coefficient of the 
55 term in (3.61) is correct. While it is not possible to be certain from these results that 
the O(S) terms are also correct, the divergence of the asymptotic results from numerical 
results arises at similar values of 5 to the symmetric case, when the 0(5) terms were 
seen to be correct. In addition, it should be noted that the next higher order term in 
(3.61) is 0(5 Q2 ), a 2 « 1.387 which is very close in magnitude to the 0(5) term retained. 
Comparison of numerical results and asymptotic solutions in the nose region suggests that 
the coefficient of the next term in the expansion (3.41) is fairly large, further suggesting 
that the 0(5 Q2 ) in (3.61) is indeed significant, even for small 5. Hence figure 3.2c and d 
shows that the asymptotic form (3.61) is at least consistent with numerical results. 

The absence of any 5 In 5 term in the expansions of the symmetric and anti-symmetric 
receptivity coefficients is somewhat surprising. However, further insight can be gained 
by considering the case arg(£) = — 7r/4 in more detail. In this special case, the operators 
AA(°),AT (3) and 77^ defined in Appendix C are real. Noting that the set of functions H(rj) 
are purely imaginary, equations (3.36) show that the functions ^,(£,77) are also wholly 
imaginary. Hence in the limit |£| — ► oo, w 3 /w 0 is imaginary and Re(c3^) = 0. That is, 
if there is a 5 In 5 term in the small-5 expansion of C s (S ), then it is wholly imaginary. 
Furthermore, we can show from the large-£ limit of w 4 /w 0 , that 

Re(4 3) ) = -|lm(r, (l) + c { 3 s) ). (3.62) 


(0) 

c a m(e/2) r ' [l + otr 1 )], 

ci a) 

4 a) -AA 1( £, 


(3.58) 


4 a) + T-r ln(£ 2 /2) - A£(2A) + f -B x -A x ln(£ 2 )), 


.(i) 


42 



These results are in agreement with the numerical results. The same results also apply 
to the anti-symmetric case. Further consideration of the large-£ limit of ic 0 and u 0 for 
arg(£) = — tt/ 4 proves that 

arg(C J (0)) = -|(l + r 1 (0 »). arg(C„(0)) = -|(i + t<°») (3.63) 

which confirms the numerics, and also explains the tt/4 difference in phase between flat- 
plate receptivity coefficients for symmetric and anti-symmetric forcing noted in chapter 
2 . 

3.6 Summary 

In §3.5, the small-5 expansions of the real and imaginary parts of the receptivity coeffi- 
cients for symmetric and anti-symmetric forcing were compared with numerical results. 
From a practical standpoint, the most important quantity is the modulus of the recep- 
tivity coefficient, since it is the variation of this that is likely to influence the position of 
the transition point. In figure 3.3, the asymptotic prediction for the modulus of the re- 
ceptivity coefficient is compared to numerical results, for symmetric and anti-symmetric 
forcing. For oblique acoustic waves ( 8 ± 0), the coefficient K a in (3.20) is often large 
compared to k 3 (see chapter 2). Moreover, as 5 — ► 0 the anti-symmetric receptivity 
coefficient C a is approximately five times larger than C 3 , the symmetric receptivity co- 
efficient. Hence C a (S ) is much more important than 67,(5 ) in determining the behavior 
of the total receptivity coefficient Ci(5). In this chapter, we have shown that C a (S ) has 
a singular structure for small 5, varying as 5a, so that small changes in this parameter 
have a significant influence on the receptivity. Moreover, the exact form of this leading- 
order behavior has been accurately calculated. The rather small range of validity of the 
asymptotic expansions at first appears somewhat disappointing. However, numerical re- 
sults show that the maximum receptivity occurs for small 5 and hence the behavior in 
this limit is of most interest for practical applications. For 5 < 0.005, excellent agreement 
is seen between asymptotic theory and numerical results, so in this parameter range the 
asymptotic expansions developed in this chapter are of quantitative as well as qualitative 
value. 

The small-5 limit is also of interest when considering the experimental results of Saric 
et al. (1995). In these wind-tunnel tests, a leading edge consisting of a machined super- 
ellipse of either 1:20 or 1:40 aspect ratio is attached to a flat plate. By considering a 
super-ellipse profile rather than an ellipse for the leading-edge section, localized receptiv- 
ity due to a discontinuity in surface curvature at the join with the flat plate (Goldstein 
1985) is eliminated. The region of dominant receptivity should then be the leading edge. 
Measurements were made for free-stream speeds 8ms -1 < U < 21ms -1 , corresponding to 
nose radius Reynolds numbers Re n = Ur n /u in the range 130 < Re n < 670 , consistent 
with the assumption Re n 1 used in our analysis. Acoustic forcing at non-dimensional 
frequencies in the range 30 x 10 -6 < e 6 < 90 x 10 -6 was considered, leading to a Strouhal 
number range 0.004 < 5 < 0.060. Our small-5 theory is clearly relevant to the lower end 
of this Strouhal number range. The results of the present analysis show that, despite the 
small values of the nose radius in these experiments, there can be significant departures 
from the flat-plate results. Many of the detailed measurements of Saric et al. are for a 
Strouhal number of 0.01. For this value of 5, symmetric disturbances lead to receptivity 
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0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 

Figure 3.3: Comparison of numerical and asymptotic results for the modulus of the 
receptivity coefficient for (a) symmetric forcing, |C,| vs. S, and ( b ) anti-symmetric 
forcing, |C 0 | vs. S. Numeric results are denoted by o, asymptotic theory is marked by 
the solid line. 
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2% greater than the flat plate value, but for anti-symmetric forcing there is a 15% re- 
duction compared to the flat-plate result. As noted earlier, despite the small numerical 
values of S, by S = 0.01 the computed value of the receptivity coefficient has already 
begun to diverge from the asymptotic results. 

For leading-edge receptivity, various measures of receptivity can be considered. We 
have chosen to define C|, the coefficient of the first Lam-Rott asymptotic eigenfunction, 
as the ‘receptivity coefficient’. Alternative definitions of receptivity level that have been 
considered in computations and experiments are values based on the Tollmien-Schlichting 
wave amplitude at the lower neutral-stability point (branch I), or on an extrapolation 
of the Tollmien-Schlichting wave amplitude back to the leading edge. Extrapolation of 
the Tollmien-Schlichting wave amplitude back to the leading edge is not advisable, since 
the Tollmien-Schlichting wave is not a valid solution of the disturbance equations near 
the leading edge. In fact, in extrapolating back to the leading edge the slowly varying 
amplitude A(£) multiplying the mode-shape function of the Tollmien-Schlichting wave 
should be considered, and A(£) becomes infinite at the leading edge, the streamwise 
velocity fluctuation behaving as £ 2n where T\ = 4- is given by (3.18). 

A receptivity level defined in terms of the Tollmien-Schlichting wave amplitude at 
the lower branch neutral-stability point does allow easy comparison with experiments 
or computations, but has the disadvantage that the results depend on the frequency 
and Reynolds number, making presentation of the results much less compact. However, 
there is a need to relate receptivity levels defined in terms of Ci to levels expressed in 
terms of branch I amplitudes. In the asymptotic theory, this is achieved by developing the 
asymptotic expression for the Tollmien-Schlichting wave in the OSE region and matching 
this expression to the first generalized Lam-Rott asymptotic eigensolution (3.18). This 
procedure is discussed by Goldstein (1983) for the flat-plate case, and the extension 
to non-zero S is discussed briefly in chapter 2. However, the asymptotic analysis is 
quite laborious, especially for the determination of A(£). An alternative approach is 
the numerical solution of the disturbance equations, using the Lam-Rott asymptotic 
eigensolution as a starting condition and continuing downstream to branch I. 

From a theoretical standpoint, the most attractive measure of receptivity is the coef- 
ficient Ci of the first Lam-Rott asymptotic eigenfunction. This Lam-Rott eigenfunction 
evolves into the unstable Tollmien-Schlichting wave farther downstream so that, when 
appropriate scaling factors related to the asymptotic matching of the LUBLE and OSE 
regions axe introduced, Ci is also the coefficient of the Tollmien-Schlichting wave. An 
extremely attractive feature of this receptivity measure is that Ci is independent of the 
physical frequency and the (asymptotically large) Reynolds number, thus providing the 
receptivity results in the simplest dimensionless form. The quantity Ci is determined by 
the unsteady flow behavior upstream of the region of instability and therefore focuses 
on the receptivity aspect of the unsteady flow development, with less influence from the 
global stability properties of the flow. Also, for fixed values of the frequency and Reynolds 
number, variations in the receptivity coefficient C\ with changes in the free-stream distur- 
bance characteristics translate directly into variations in the Tollmien-Schlichting wave 
amplitude. 
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Chapter 4 


Receptivity for a leading edge with 
aerodynamic loading 

4.1 Introduction 

The transition of a boundary-layer flow from laminar to turbulent is strongly influenced 
by free-stream disturbances. The process by which free-stream disturbances generate 
instability waves in the boundary layer is known as receptivity (Morkovin 1969). Energy 
is transferred from long-wavelength disturbances in the free stream to the much shorter 
wavelength Tollmien-Schlichting waves in the boundary layer. This transfer requires a 
wavelength conversion mechanism, which is usually produced by non-parallel mean flow 
effects due to short-scale streamwise variations in the mean flow. These streamwise gradi- 
ents may arise from slow, viscous boundary-layer growth associated with the leading-edge 
region (Goldstein 1983, Hammerton & Kerschen 1996, 1997), or from localized regions 
farther downstream in the boundary layer, such as regions with sudden changes in sur- 
face geometry (Goldstein 1985, Choudhari and Kerschen 1990) or marginally separated 
regions (Goldstein, Leib & Cowley 1987). Additional background information on recep- 
tivity can be found in Goldstein & Hultgren (1989) and Kerschen (1990). 

Leading-edge receptivity was first considered for the Blasius boundary layer on a flat 
plate (Goldstein 1983). A high Reynolds number asymptotic analysis was formulated for 
an incompressible, two-dimensional flow of free-stream speed Uoo, together with a small- 
amplitude time-harmonic perturbation of frequency u>, equivalent to the incompressible 
limit of an acoustic wave propagating parallel to the mean flow. Two streamwise regions 
enter the analysis, one region of 0(Uoo/uj), where the inviscid pressure field and slip ve- 
locity induced by the free-stream disturbance drives the unsteady motion in the bound- 
ary layer, and a second region farther downstream at distance O^uuj/U ^) - 3 U^/uj), 
where the disturbance is governed by the triple-deck structure, corresponding to the high 
Reynolds number asymptotic form of the Orr-Sommerfeld equation (OSE) in the vicinity 
of the lower branch. 

In the first region of Goldstein’s analysis, assuming that the amplitude of the unsteady 
disturbances in the free-stream is small compared with the mean flow, the unsteady 
disturbance in the boundary layer is governed by the linearized unsteady boundary layer 
equation (LUBLE). Far downstream in the LUBLE region, the solution consists of a 
Stokes wave, and a set of asymptotic eigenfunctions that contain velocity but not pressure 
fluctuations. These asymptotic eigenfunctions are equivalent to those obtained by Lam & 
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Rott (1960, 1993) and Ackerberg & Phillips (1972). The asymptotic analysis for distances 
far downstream determines the form of the eigenfunctions, but not their coefficients C*. 

The wavelengths of the asymptotic eigenfunctions shorten progressively with distance 
downstream. Eventually, the self-induced pressure field associated with the displacement 
thickness of each asymptotic eigenfunction becomes significant, and the triple-deck struc- 
ture replaces the LUBLE as the correct asymptotic approximation to the Navier-Stokes 
equation. The first asymptotic eigenfunction of the LUBLE matches on to the Tollmien- 
Schlichting wave solution of this triple-deck region. Thus, the form of the free-stream 
disturbance and the geometry close to the nose influence the amplitude of the Tollmien- 
Schlichting wave only through the coefficient C\ of the first asymptotic eigenfunction. 
For this reason, C\ is known as the ‘receptivity coefficient’. The numerical value of the 
receptivity coefficient can not be determined by asymptotic methods. Instead it must be 
extracted from numerical solutions of the LUBLE, by comparison with the asymptotic 
form far downstream in the LUBLE. This was accomplished for the flat-plate case by 
Goldstein, Sockol & Sanz (1983) and Heinrich & Kerschen (1989). 

The analysis of Goldstein for a flat plate was the first theoretical description of the 
fundamental mechanisms of leading-edge receptivity. However, aerodynamic bodies de- 
signed for subsonic flow generally have finite thickness distributions with a parabolic lead- 
ing edge. Hammerton & Kerschen (1996, 1997) considered a thin, symmetric airfoil at 
zero angle of attack with a plane acoustic wave incident at arbitrary angle, and examined 
the influence of the leading-edge geometry of the body on receptivity. The structure of 
the development of the instability remains similar to that for a flat-plate; the asymptotic 
eigenfunctions of the LUBLE are modifications of the Lam-Rott eigenfunctions which 
take account of the effects of mean pressure gradient and surface curvature. Attention 
was focused on the variation of the receptivity level with the nose radius of the body 
and the incidence angle of the acoustic field. The radius of curvature of the leading-edge 
of the airfoil, r n , enters the analysis as a Strouhal number, S = ujr n /U !X> . The results 
revealed that, for disturbances parallel to the mean flow, a small increase in the recep- 
tivity for very small S is followed by a rapid decrease in the level of receptivity so that, 
when S = 0.3, the receptivity is reduced to approximately 15% of the flat-plate value. In 
addition, for unsteady disturbances propagating at an angle to the mean flow, the over- 
all level of receptivity is dominated by the response to the component of unsteady slip 
velocity which is anti-symmetric about the leading-edge. Inclusion of thickness effects 
goes some way towards modeling experimental investigations (Saric, Wei, Rasmussen & 
Krutekoff 1995; Saric & White 1998) and allowing comparisons with related numerical 
simulations (Reed & Lyttle 1998). In wind tunnel tests reflections of acoustic waves from 
the side walls can lead to significant anti-symmetric flow about the leading edge (see 
Kerschen 1990 and Saric & White 1998). 

In this chapter the receptivity of a cambered airfoil at an angle of attack to a mean 
flow is considered. The receptivity problem for this physical situation is considerably 
more intricate than the examples analyzed up to this time. First, there is receptivity 
associated with the mean stagnation point on the lower surface, as well as receptivity due 
to the strongly non-parallel flow around the leading edge. Moreover, at non-zero angles 
of attack a minimum appears in the steady wall shear on the upper surface. As the angle 
of attack of the airfoil is increased, the value of the minimum wall shear decreases until 
at some critical angle the minimum wall shear reaches zero, which for a two-dimensional 
flow corresponds to flow separation. Hence a regime exists where the angle of attack 
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is close to the critical angle and there is a region of marginal separation on the upper 
surface of the airfoil. Here the receptivity mechanism analyzed by Goldstein et al (1987) 
is relevant. 

The structure of the chapter is as follows. In §4.2 a high Reynolds number asymptotic 
analysis (e 6 = vu>/Ul\ ^ ^ 1) is formulated for an incompressible, two-dimensional flow 
about a thin, cambered airfoil. In §4.2.1, the steady slip velocity on the airfoil surface is 
determined using thin-airfoil theory, followed in §4.2.2 by the derivation of the equation 
governing the unsteady boundary-layer flow in terms of suitable coordinates. Assuming 
that the unsteady perturbations in the free-stream are small, and letting e — ► 0, the 
unsteady disturbances in the boundary layer are governed by the LUBLE. 

The structure of the solution in the LUBLE region is analyzed in §4.3. When consid- 
ering the subsequent development of instabilities in the boundary-layer flow, the asymp- 
totic form of the unsteady solution far downstream in the LUBLE region is required. 
The derivation of the asymptotic form of both the steady and unsteady components of 
the solution is a necessary part of the presentation, but details are given in Appendices 
D and E instead of the body of the chapter in order to avoid interrupting the description 
of the evolution of the boundary-layer disturbance. 

In §4.4, the inviscid pressure field and slip velocity produced by the interaction of a 
free-stream acoustic wave with an airfoil are determined, in order to allow calculation 
of the interaction between the free-stream and the boundary layer. The exact form 
of the unsteady slip velocity in the vicinity of the leading edge is determined by the 
global solution about the airfoil. This depends on the magnitude of the reduced acoustic 
frequency k = ub/c, where b is the airfoil semi-chord and c is the speed of sound. Here we 
present results for the limiting cases k •C 1 and k 1, when relatively simple expressions 
for the slip velocity can be obtained. The definition of the receptivity coefficient for 
this particular geometry is also provided. Finally in §4.5, numerical solutions of the 
steady boundary layer equation and the LUBLE are obtained, and compared with the 
asymptotic eigenfunctions of §4.3. This allows determination of the receptivity coefficient 
as a function of leading-edge geometry, airfoil shape and angle of attack, and of the 
characteristics of the free-stream acoustic wave. 


4.2 Formulation 

We consider a thin, cambered airfoil of chord 2b at an angle of attack q to a uniform 
flow of speed U^. A plane acoustic wave of frequency u>, propagating at an angle 6, is 
assumed to be incident on the airfoil as illustrated in figure 4.1. Two-dimensional, low 
Mach number flow is considered. Since the Mach number is small, the mean flow can 
be analyzed using incompressible theory. For the unsteady component of the flow, most 
features of interest can also be analyzed using incompressible theory. The influence of 
compressibility on the unsteady component of the flow is discussed in §4.4. The Reynolds 
number is assumed large, so the flow field is inviscid and irrotational everywhere except 
in the vicinity of the airfoil surface. The behavior of the flow in the boundary layer 
adjacent to the airfoil surface depends crucially on the pressure gradient imposed by the 
outer inviscid flow. Therefore, relevant features of the inviscid mean flow are presented 
in §4.2.1. The equations governing the boundary-layer flow are formulated in §4.2.2. 
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Figure 4.1: An illustration of the physical situation of interest: a thin, cambered airfoil of 
chord 26 is at angle 5 to a uniform flow of speed U 0 0 , with a plane acoustic wave incident 
at an angle 6. 


4.2.1 Inviscid Mean Flow 

We introduce Cartesian coordinates (x, y) normalized by the airfoil semi-chord 6, with 
the origin located at the airfoil leading edge and the x- and y-coordinate axes parallel 
and normal to the airfoil chord line at the leading edge, as illustrated in figure 4.1. This 
particular choice of coordinate system proves convenient for the present analysis, in which 
attention is focused on the leading-edge region. The camber and thickness of the airfoil 
are assumed of 0(5) where 5 1. Explicitly factoring out the scaling parameter 5, the 

airfoil surface is defined by 

y = 5(-n(x) ± s(x)), 0 < x < 2, (4.1) 

where the ± sign applies on the upper and lower surface, respectively, and the camber 
and thickness functions n(x) and s(a:) are defined in figure 4.1. The angle of attack 5, 
measured with respect to the axes defined above, is also assumed to be small (to avoid 
separation) and to be comparable to the thickness and camber; hence we write a = 6a. 
The inviscid mean flow past the airfoil can then be calculated using thin-airfoil theory. 
Specifically, the no-penetration boundary condition on the upper and lower surfaces of 
the airfoil can be linearized and transferred to y = 0*, respectively. The mean-flow 
perturbation due to the presence of the airfoil is found to have the y-components 


6/ 00 6 (— a + n'(x) ± s'(a:)) on 0 < x < 2, y = 0*, (4.2) 

where n'(x) and s'(x) denote the derivatives of these functions with respect to x. The 
mean slip velocity on the airfoil surface, C/ 3 , immediately follows (Cheng & Rott 1954), 


a'fri) . 

2:1—2: 


2 — xl 2 f 2 [ 2:1 ]5n'(2: 1 ) — a 


2 - x, 


Xi — x 


dx ! + 0(6 2 ). (4.3) 
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A circulatory component has been included in (4.3), its magnitude fixed by the Kutta 
condition at the trailing edge. 

The thin-airfoil expansion (4.3) is invalid in the vicinity of the leading edge, and must 
be supplemented by a local expansion. For an airfoil with a rounded leading edge of 
(dimensional) radius r„, the small argument expansion of the thickness distribution has 
the form 

s (x) = SiX l/2 + s 2 x + s 3 x 3/2 + 0(x 2 ) (4.4) 

where Si = (2r n /6 2 b ) 1 ^ 2 . The coefficient of the leading term in (4.4) must be 0(1), 
implying that r n = 0(S 2 b). The airfoil nose radius r n is the natural length scale in the 
leading-edge region. We therefore introduce a parabolic coordinate system (<f, fj) based 
on this length scale. The relationship between the (x, y) and (£, fj) coordinates is 

x+iy = \jl^ +if i ) 2 + 1 ]- ( 4 - 5 ) 


In terms of these new coordinates, the airfoil surface is given by fj = 1 + |<5s 2 £ 3 /(l + £ 2 ) + 
0(6 2 ) . Thus, at leading order in the local coordinates, the airfoil surface is defined by the 
parabola fj = 1. The inviscid flow in the vicinity of the parabolic nose can be obtained in 
the form of a complex potential by means of conformal mappings, giving a slip velocity 

£ 4(0 - U t ^| 2 + ^1/2 + + 1)1/2) • ( 4 - 6 ) 

The first term of (4.6) corresponds to symmetric flow past the nose, while the second 
term is an antisymmetric flow around the nose, from the lower surface to the upper, 
when the parameter y is positive. The parameters U e and y in (4.6) are determined by 
matching with the thin-airfoil expression (4.3). 

In order to match (4.3) and (4.6), the small-rr behavior of (4.3) and the large-f 
behavior of (4.6) must be determined. For i«l, (4.3) takes the form 


U.(x) 


= l + — ( ±na e \/2x 1/2 — s 2 lnx + -^= + s 2 In 2 
n \ v2 


- I* [— _ 

Jo 


[cbi 2x\ /2 


ckci 

Xi 


(4.7) 


where a e , the effective angle of attack for the leading-edge region taking into account the 
camber of the entire airfoil, is given by 


1 c 

oc e — ct / 

7 r Jo 


n'(x) 

M 2-*,)] 1 / 2 


dar!. 


(4.8) 


For f » 1, (4.6) takes the form U s (£) = U e + U e y/£ + 0(l/£ 2 ). The matching of the 
symmetric components of these expressions is identical to that discussed in chapter 2. 
For leading-edge shapes that contain a wedge component (s 2 7^ 0), the approximation 
(4.6) is valid only at 0(1) in the thickness parameter S and the matching then shows 
that U e = U oo. For leading-edge shapes, such as that of a Joukowski airfoil, where the 
even coefficients in (4.4) vanish, the leading-edge region is also parabolic at 0(<5), and 
the matching then gives 



(4.9) 
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The 0(6) term in (4.9) is essentially a correction to the ‘free-stream speed’ in the local 
leading-edge region, due to the flow perturbation created by the thickness distribution 
of the entire airfoil. Next, matching of the antisymmetric components of the local and 
thin-airfoil expansions then gives the angle-of-attack parameter 


^ = 2Qe V^T’ (4 ' 10) 

which is 0(1) since r n = 0(5 2 b). The stagnation point on the lower surface is given by 
f = ~n, which corresponds to x = 2 (6a e ) 2 , a familiar result for thin bodies. 


4.2.2 Boundary-layer flow 

The boundary-layer flow in the vicinity of the leading edge is also analyzed most conve- 
niently in parabolic coordinates. For the mean boundary layer, the nose radius r n remains 
the appropriate length scale in the streamwise direction. However, for the unsteady flow, 
a second streamwise length scale, UJ w, is also relevant. Specifically, Goldstein’s (1983) 
analysis of leading-edge receptivity for the flat-plate boundary layer showed that the 
receptivity process takes place in a region where the distance from the leading edge is 
0(U e /u). To focus on this region, we introduce new coordinates, 

£ = V = S*fj, (4.11) 

where 

Q w 1 Tl 

~~ul 

is a Strouhal number based on the airfoil nose radius. 

Since the Reynolds number is assumed large, viscosity is important only in 
boundary layer adjacent to the body surface. To analyze the boundary-layer flow, 

i = f > f) - S* = e 3 r), 

where the small parameter 

‘ = jji « 1 (4.14) 

is the reciprocal of the Reynolds number based on the disturbance length scale U e /ui. 
The analysis presented here considers the small-e limit but with S = 0(1). The non- 
dimensional vorticity equation, expressed in terms of the streamfunction 'I' (which has 
been normalized by e 3 U 2 /u), then becomes 


(4.12) 

a thin 
we set 

(4.13) 






d(Z,r)) 


mjn _ f 6 


h 2 




h P 


■ + 

r^i 

[H 2 \ 

+6 12 

r*« 


1 

H 2 






d^/H 2 , £) 

d(^T]) 



(4.15) 


where H = (^ 2 + 5)T The streamfunction 'F contains both the mean and time-dependent 
components of the flow. Equation (4.15) is exact, except for approximation of the metric 
coefficient (£ 2 + fj 2 ) by H . The viscous flow satisfies the no-slip boundary conditions on 
the body surface, 

'F = 'F„ = 0, on ?7 = 0. (4.16) 
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At the outer edge of the boundary layer, the viscous flow matches to the inviscid slip 
velocity (non-dimensionalized by U e ), 

-* U 3 (£,t), as (4.17) 

For 0(1) values of £, the terms on the right-hand side of (4.15) can be neglected, 
leading to the unsteady boundary-layer equation. The unsteady component of the flow, 
a small perturbation to the mean flow, then satisfies the linearized unsteady boundary- 
layer equation (LUBLE). It is shown in §4.3 that the solution of the LUBLE contains 
components whose wavelengths progressively shorten with distance downstream. Thus, 
for the unsteady component of the flow, terms on the right-hand side of (4.15) become 
significant when £ = 0(e -1 ), and the correct asymptotic approximation to (4.15) then 
takes on the triple-deck structure. The asymptotic matching of these two streamw'ise 
regions is discussed briefly at the end of §4.3. 

4.3 Analysis of receptivity region 

In this section, we consider the region of the boundary layer where the receptivity takes 
place. The streamwise length scale for this region is U e /u>. Thus, the development of 
the viscous flow in this region is governed by the form of (4.15) corresponding to the 
limit e — » 0, £ = 0(1). Subsequently, the behavior of this solution at large values of 
£ will be considered, in anticipation of matching with the Orr-Sommerfeld region that 
exists farther downstream in the boundary layer. For £ = 0(1) and e — > 0, terms on the 
right-hand side of (4.15) can be ignored. Integrating once with respect to 77 then gives 
the unsteady boundary-layer equation, 

H 2 ^ + (¥„% - ¥„,¥*) - £/r 2 tf„ 2 - ^ = - H 2 P t j, (4.18) 

where the pressure gradient F((£, t ) is obtained by matching to the outer (inviscid) flow, 

(4.19) 

Up to this point our analysis has paralleled that of chapter 2, utilizing nearly the 
same notation. For the symmetric mean flow case considered by chapter 2, the unsteady 
boundary-layer region was analyzed using the streamwise coordinate £ (which is based 
on the length scale U e /u>). This allowed easy comparison with the flat-plate analysis 
of Goldstein (1983). However, for the non- zero angle-of-attack case considered here, no 
corresponding flat-plate analysis is possible because the boundary layer would separate 
at the leading edge. In addition, the presence of airfoil thickness, camber and angle- 
of-attack significantly complicates the algebra. The analysis of the mean flow is most 
naturally carried out in terms of a streamwise coordinate based on the nose radius r„. 
Noting the parabolic nature of the governing equations, one can see that the mean-flow 
stagnation point is the natural origin for the streamwise coordinate. Thus we introduce 
the new streamwise coordinate, 


w = £ + n = S *£ + /x, 


(4.20) 
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where // is the angle-of-attack parameter defined in (4.10). The mean-flow stagnation 
point is now given by w = 0, the nose corresponds to w = n and w is scaled on the nose 
radius r n rather than the disturbance length scale UJ u. 

Because the mean flow is expressed most simply in terms of the streamwise coordinate 
w, it turns out that w is also the most convenient coordinate for analysis of the unsteady 
flow. Our analysis assumes S = 0(1), so that w = 0(1) also corresponds to f = 
0(1). Later in this section we analyze the behavior of the solution for large £, in order 
to examine the development of the Lam-Rott asymptotic eigenfunctions and the 
eventual evolution of the first eigenfunction xpi into the Tollmien-Schlichting wave when 
£ = 0(l/e). This large-£ analysis will be presented in terms of large w for algebraic 
simplicity, but we emphasize that it is the quantity S l ^ 2 w which is required to be large. 

For typical values of the camber distribution and angle of attack, the mean aero- 
dynamic loading parameter n is usually positive. VVe examine the influence of mean 
aerodynamic loading on the receptivity for both surfaces of the airfoil. However, rather 
than considering the cases w > 0 (upper surface) and w < 0 (lower surface) separately, 
the analysis can be presented more concisely by considering w to be positive and taking 
/x > 0 to give upper surface results and [i < 0 for lower surface results. 

For a cambered airfoil at non-zero angle-of-attack, the slip velocity is given in terms 
of the new coordinate w by 


~ w 

U a (w, t) = — + A u 3 e~'\ (4-21) 

where h 2 = 1 4- (to — /x) 2 . Here the steady contribution U 3 — w/h follows from (4.6), while 
the time-dependent component u 3 (w) depends on the particular form of the free-stream 
disturbance, as discussed in §4.4. Since we are concerned only with small-amplitude free- 
stream disturbances, that is A <t; 1, the steady and unsteady components of the flow 
field can be analyzed separately. Thus, the streamfunction within the boundary layer 
can be written in the corresponding form 

^(it;, 7]) = S^w(f}(w, tj) + Aip(w, T))e~' 1 , (4.22) 

where the factor S^w has been extracted from the mean component of the streamfunction 
in order to simplify the matching condition as t] — » oo. The steady boundary-layer flow 
then satisfies 

$T)TIT] 4" 07^)0 "b ^(0TP70U) 0T/077ti>) /^(0IJ 1) = 0, (4.23) 

where the pressure gradient parameter j3(w\ ji) is given by 


0 = 


1 — fiw + n 2 

h? 


(4.24) 


together with boundary conditions 

0 = 0,, = 0 at 77 = 0 and 4> n —* 1 as 77 — > oo. (4-25) 

The time-dependent component of the boundary-layer flow satisfies the linearized 
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unsteady boundary-layer equation (LUBLE), 


JF( 0 ) = SVi 2 


dp 
d w' 


?{4>) = trm + (0 + wtwWrm + iSh 2 + 


2 , W 2 - /X 2 - 1 


h? 


<Pi 7 W(prj W Ipjrj 


~\~'W{(p T jrj'lpUJ (prj^Prjw ) 


d P ! cu 1 - K w ~ A0\ wdu s 

^={' sh P 

with boundary conditions 

ip = ipT) = 0 at 77 = 0 and ip v — > S%hu s as 77 — ► 00 . 


(4.26) 


(4.27) 


These equations must be solved numerically for each value of S and /x, and for each 
different free-stream disturbance, although various simplifications arise as will be seen 
later. The numerical solutions are described in §4.5. However, as discussed in §4.1, we 
are primarily interested in the solution of these equations in the large-£ (or large- 7 / 7 ) limit, 
where a component of the solution to the LUBLE matches onto the Tollmien-Schlichting 
wave solution of the Orr-Sommerfeld equation. In the next two sub-sections, we develop 
large-in asymptotic expansions for 0 ( 117 , 77 ) and 0 ( 1 x 7 , 77 ). 


4.3.1 Large w expansion of the steady boundary-layer equation 


In order to determine the form of the asymptotic eigensolution of the LUBLE that 
matches onto the TS wave far downstream, it is found that the expansion of the mean 
flow must be determined up to, and including, terms of 0(w~ 3 ). The large-tu expansion 
of /? is 


W W z 7X7 J \W 4 / 


(4.28) 


For /x > 0 an adverse pressure gradient arises on the upper surface (w > 0); this adverse 
pressure gradient decays to zero much more slowly than the favorable pressure gradient 
arising in the symmetric mean-flow case (/x = 0). Guided by the p = 0 case (Van Dyke 
1964), we find that the asymptotic expansion of the mean flow far downstream is given 
by 


0 ( 117 , 77 ) ~ F( 77 ) 


miv) | Mrj) logw 2 Psjrj) 
w w 2 w 2 


f^Mv) log*" 2 h f^o(jj) +0 f l \ 

W 3 7X7 3 \W~ r2 ) 


(4.29) 

where the fractional power 72 = 3.774 arises as an eigensolution (Libby & Fox 1963). 
Here F(rj) is the Blasius solution and the pi(r)) are determined by the set of equations 
£»(p») = di, given in Appendix D. The parameter /x appears in the p,. As for the 
symmetric mean-flow case described in chapter 2, the large-7X7 asymptotic form of the 
mean flow is not fully determined by the local conditions far downstream. A coefficient 
B i(/x) enters the 0(xx7 -2 ) and 0(w~ 3 ) terms, the value of which is dependent on conditions 
close to the leading edge. Hence £?i(p) must be determined by comparing the asymptotic 
form (4.29) with the numerical solution obtained by integrating (4.23) from the stagnation 
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point w — 0. The asymptotic form (4.29) for the mean flow is applicable to both the 
surfaces of the airfoil; typically fj, > 0 corresponds to the upper surface and n < 0 
corresponds to the lower surface. 

In order to determine the large-u; solution for the unsteady component of the flow, 
the behavior of the mean flow close to the wall is required. From the results of Appendix 
D, this is given by 


where 


4> 


rfU^w) | rfU^w) 


rfFp 

5! 


//2 


+ 0 ( 


w' n ’ w 


(4.30) 


uiM 
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(4.31) 

The functions Uq(w ) and Uq(w) are the asymptotic expansions of the shear stress and 
the curvature of the mean boundary-layer profile at the wall, accurate to 0(w~ 3 ). The 
numerical values of the constants Fq, Ai, ho, j 0 and ko are given in Appendix D. 
Numerical results for B i(/z), obtained by comparing the numerical solution of (4.23) 
with the asymptotic form (4.29), are presented in §4.5. 


4.3.2 Large S l ^ 2 w expansion of the linearized unsteady 
boundary- layer equation 

We next consider the evolution of the unsteady component of the flow, which is governed 
by (4.26). As in the symmetric mean-flow case considered in chapter 2, far downstream 
in the LUBLE region [S l ^ 2 w 1) the unsteady component of the flow consists of 
a particular solution, ib v , determined entirely by the local conditions far downstream, 
together with a set of asymptotic eigenfunctions, 


Ip{w, 77; S, n) = ip p (w, 77; S, n) + Ci(S, fj)ipi(w, 77; S, n). (4.32) 

i 


The particular solution is a generalization of the classical Stokes layer solution, driven by 
the local value of the unsteady pressure gradient (4.26). The asymptotic eigenfunctions 
ipi are generalizations of the Lam-Rott eigenfunctions (Lam &: Rott 1960), taking account 
of the non-Blasius mean flow. A second set of asymptotic eigenfunctions for the LUBLE 
were derived by Brown and Stewartson, but their relationship to Tollmien-Schlichting 
wave development has not been established. The eigenfunctions are ‘asymptotic’ because 
they exist only for S l ^ 2 w 1, where the mean flow takes on a slowly varying character. 
The tpi depend on the mean-flow characteristics far downstream, but are independent of 
the unsteady free-stream disturbance in the downstream region. The coefficients C* of 
the asymptotic eigenfunctions are determined by the characteristics of the unsteady free- 
stream disturbance and the mean boundary-layer flow in the region nearer the leading 
edge where S x ^ 2 w = 0(1). One of the asymptotic eigenfunctions, which we label ipi, 
matches on to the growing Tollmien-Schlichting wave in the Orr-Sommerfeld region 
farther downstream, where S l ^ 2 w = 0(l/e). It is only through the receptivity coefficient 


55 



Ci that the unsteady free-stream disturbance influences the amplitude of the Tollmien- 
Schlichting wave. In the remainder of this chapter, most attention is focused on this 
particular eigenfunction and on its coefficient C 

The form of the general asymptotic eigenfunction ipi is obtained in Appendix E. The 
asymptotic eigenfunctions have a two-layer structure, with a new inner layer of width 
T] = 0(S~ l / 2 w~ l ). The leading-order expression for the general eigenfunction ipi is given 

by 


F^pfe ( Sw 2 \ T ' 


exp (T} 0 \w))f^’(n), n = Siwtj = 0(1) 


r(0) 


W = 


/“ Ai(z)dz V 2 
■wS* ( W ^ 


(4.33) 


exp (T> >) F\r ) ), 


T} = 0( 1 ) 


where pi is the z-th root of Ai'(— p) = 0, and T[°\w) and fj°\n) are given by (E.20), 
(E.21) and (E.30), respectively. 

In the next section, numerical solutions for the LUBLE are compared to the asymp- 
totic eigenfunction in order to extract values of the receptivity coefficient C\. The 
two most convenient points of comparison are the unsteady component of the wall shear, 
= 0), and the oscillating boundary-layer thickness which is related to the value of 
the streamfunction far from the wall, xp(r} — > oo). For the first eigensolution, p x = 1.01879 
and (4.33) gives 

/ c. 2\ T i + 1 

= 0) ~ 0.4356(1 + i) I — ) exp^^u;)) [l +0 (nT 72 * 3 )] , (4.34) 


ipi(r) — > oo) ~ wS 2 exp(Ti 0) (iu)) |^1 + O (w 72+3 )] , (4.35) 

with 

T\ = -0.6921 + 1.9878(1 - 6.182/x 2 )i 5 (4.36) 

The exponent for the first neglected term in (4.34) and (4.35) is —72 + 3 = —0.774. 
This term is followed by a sequence of closely-spaced correction terms of 0(w~ l In 2 iu), 
0(w _1 lnw) and O(io -1 ), as can be seen from the analysis presented in Appendix E. 

Before proceeding to numerical solutions of (4.26), we first consider the validity of 
the LUBLE. While the eigenfunctions obtained here are uniformly valid solutions for the 
LUBLE as f — » 00, they are not uniformly valid large-£ solutions of the full equations gov- 
erning the development of the viscous flow. Since T oc £ 3 at leading order, the wavelengths 
of the eigenfunctions decrease with distance downstream, increasing the importance of 
terms involving streamwise derivatives. When £ = 0(l/e), terms on the right-hand side 
of (4.15), which were neglected in forming the LUBLE, now become significant. An ir- 
rotational layer outside the mean boundary layer, driven by the oscillating displacement 
thickness of the eigenfunction, must then be considered, and the pressure gradient im- 
posed by the motion in this outer layer appears in the leading-order equations governing 
the wall layer. This coupled viscous-inviscid interaction has the triple-deck structure, 
corresponding to the small-e asymptotic approximation to the Orr-Sommerfeld equation 
in the vicinity of the lower branch. A complete treatment of the linear development of 
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the instability wave would require an asymptotic solution for the Tollmien-Schlichting 
wave of the Orr-Sommerfeld equation, taking account of the mean pressure-gradient dis- 
tribution on the airfoil. However, for 0(1) values of 5, with e and 6 small quantities of 
the same order, the Orr-Sommerfeld region is influenced by the full airfoil camber and 
thickness distributions, necessitating a general development that does not seem justified 
in the present context. 

In this chapter we restrict ourselves to examining only the receptivity process, and 
the remainder of the chapter is concerned with determining the receptivity coefficient C\ 
as a function of the nose radius parameter S, the aerodynamic loading parameter /z, and 
the characteristics of the free-stream disturbance. This is accomplished by comparing 
numerical solutions with the asymptotic solutions obtained above. 


4.4 Free-stream Disturbances 

In order to calculate the unsteady flow in the boundary layer, the inviscid flow field which 
drives the unsteady motion in the boundary layer must be determined. In this section we 
consider the inviscid flow field produced by the interaction of a free-stream acoustic wave 
with the airfoil. The acoustic wave is assumed to be incident on the airfoil at an angle 6 
with respect to the airfoil chord, as illustrated in figure 1. The slip velocity and surface 
pressure fields generated by this interaction drive the unsteady motion in the boundary 
layer, leading to the generation of a Tollmien-Schlichting wave. Since the outer inviscid 
flow is irrotational, the unsteady pressure field is easily related to the unsteady velocity 
field. Thus, we present results for the unsteady slip velocity in the leading-edge region. 

For a low Mach number flow, the acoustic wavelength 27 xc/u is long compared to the 
hydrodynamic length scale U e /u>. Thus, outside the boundary layer, the unsteady flow in 
the vicinity of the leading edge is incompressible and irrotational. Potential flow theory 
then shows that this local flow has the form 

+ Ka(0)S-^ (4.37) 

where h 2 = 1 + (w — /z) 2 . Noting that w — fi = f , the parabolic coordinate with origin 
at the leading edge, it is seen that ( w — /z)//z and 1/h correspond to purely symmetric 
and anti-symmetric flow about the leading edge, respectively. The coefficients k 3 and K a 
multiplying these eigenfunctions are independent of the nose geometry, but depend on 
the free-stream disturbance, being determined by global features of the unsteady flow. 
The S~ 2 factor is included in the anti-symmetric term in order to retain consistency of 
notation with chapter 2. Expressions for the coefficients k 3 and n a can then be taken 
directly from chapter 2. 

The nature of the unsteady interaction of the acoustic wave with the airfoil depends 
upon the magnitude of the acoustic reduced frequency k = ub/c, where c is the speed 
of sound in the undisturbed medium. For extremely low Mach numbers, the acoustic 
wavelength is long not only compared to the hydrodynamic length scale, U e /ui, but also 
compared to the airfoil chord, i.e. fc -C 1. The unsteady interaction of the acoustic 
wave with the airfoil can then be analyzed using the classical unsteady airfoil theory for 
incompressible flow, giving 

k 3 = cos 0, K a =a* sin# (4.38) 
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where a — u>b/U e is the aerodynamic reduced frequency. If we consider the other limiting 
case, ub/c ;» 1, the acoustic wavelength is then short compared to the airfoil chord, and 
the interaction of the acoustic wave with the leading edge can be analyzed by taking 
the airfoil chord to be semi-infinite. As in the first case, on the scale of the acoustic 
wavelength, the airfoil appears at leading order as a zero-t hick ness plate and the problem 
reduces to the classical Sommerfeld diffraction problem. Matching to the local solution 
in the vicinity of the leading edge, then gives 

k s = cos 0, K a = -sin^0, (4.39) 

(n M e )a 


where M e — U e /c is the Mach number of the mean flow, which has been assumed small. 

As described in chapter 2, writing ip = K 3 ip s + K a ip a , and substituting into (4.26), it 
follows that the receptivity coefficient for the upper surface is given by 

Cx = K,(e) C,(S,fi) + K a (d) C a (S,n), (4.40) 


where C 3 and C a are extracted from the solutions of 




?(tpa) = i Sh 2 - 


1 + /r - w 2 
h? 


(4.41) 


respectively, with T defined in (4.26). Here C 3 is a measure of the receptivity due to that 
component of the free-stream disturbance which leads to symmetric flow in the vicinity 
of the leading edge, and C a is a measure of the receptivity arising from that component of 
the free-stream disturbance which leads to anti-symmetric flow about the leading edge. 

The separation into components of the free-stream disturbance which are symmetric 
and anti-symmetric about the leading edge is a natural choice for consideration of airfoil 
thickness effects in the absence of aerodynamic loading, since in this case the symmetric 
and anti-symmetric components exhibit quite different dependence on the nose radius 
parameter S = ur n /U e . Specifically, when (jl = 0, the symmetric component of the slip 
velocity (4.37) varies monotonically from zero at£ = 0tolas£ — ♦ oo, for all values of 
S. In contrast, for the anti-symmetric component of (4.37), the magnitude of the slip 
velocity near the leading edge is proportional to S -1 / 2 and hence becomes large when 
the nose radius is small, due to flow around the ‘sharp’ edge. The mean flow stagnation 
point is at £ = 0 when /z = 0, so that this difference in behavior is also reflected in the 
initial conditions for the numerical solutions of (4.41) for ip 3 and ip a . The behavior of the 
corresponding receptivity coefficients, C 3 (S, 0) and C a (S, 0), was discussed in chapter 2. 

When aerodynamic loading and airfoil thickness effects are combined, the situation 
is more complex. The mean flow stagnation point is then at £ = — n (or w = 0), so that 
both components of (4.37) have a finite unsteady slip velocity across the stagnation point 
where the solutions of (4.41) are started. In this case, an alternative separation of the 
unsteady slip velocity into components that are locally symmetric and anti-symmetric 
about the mean flow stagnation point can be considered. This takes the form 


w , . 1 

= K a — + («a - HS*K 3 ) 


S2h 


(4.42) 
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where the first component vanishes at the mean flow stagnation point while the second 
component is finite. The corresponding decomposition of the receptivity coefficient is 

Ci = k,C 3 sp + (k q — vSiK a )C aSP , (4.43) 

where C 3 sp and C a sp are extracted from the numerical solutions of 

Of course, the symmetry and anti-symmetry about the mean flow stagnation point of 
the two components of the slip velocity (4.42) is only a local property, since h is not a 
symmetric function of w. However, for large values of S, the region of receptivity should 
be concentrated near the stagnation point, so that one might expect the differences in the 
local characteristics of these two components of slip velocity to be reflected in differences 
in the behavior of the corresponding receptivity coefficients, C 3 sp and C a sp. This is 
explored in the next section. 

The two sets of Receptivity Coefficients are related by 

C 3 = C 3 SP — {JiS^CaSP, C a = CaSP ■ (4.45) 

In the numerical computations, it proves convenient to solve (4.44)for ip 3 sp and ip a SP, 
thereby obtaining C 3 sp and C a sp, and to calculate C 3 and C a from (4.45). In the next 
section, the variation of the receptivity coefficients with the angle-of-attack parameter 
and the nose-radius parameter S is investigated numerically. 



4.5 Numerical Results 

In §4.3, generalizations of the Lam-Rott asymptotic eigenfunctions were obtained, taking 
account of the pressure gradient due to leading-edge thickness and mean aerodynamic 
loading. The asymptotic analysis determines the form of these eigenfunctions, but not 
their coefficients Ci. The coefficients are determined by the characteristics of the free- 
stream disturbance in the LUBLE region, f = 0(1). The first of these eigenfunctions 
is the precursor of the Tollmien-Schlichting wave. Thus, the free-stream disturbances 
influence the amplitude of the Tollmien-Schlichting wave only through the coefficient 
Ci. It appears that the receptivity coefficient C\ for a particular free-stream disturbance 
can be determined only by numerically solving the LUBLE over the full range of £ and 
examining the behavior for large £. 

Numerical solutions are required for both the mean flow, governed by the nonlin- 
ear partial differential equation (4.23), and the linearized disturbance, governed by the 
LUBLE (4.26). For convenience, the numerical solutions are computed in terms of the 
variable w rather than £ (see (4.20)). Since the governing equations are parabolic, the 
solutions were computed by marching away from the mean flow stagnation point (w — 0) 
using a Keller Box scheme (Keller & Cebeci 1970). The initial conditions at w = 0 cor- 
respond to the steady Hiemenz flow and its quasi-steady linear perturbation. As noted 
earlier, in order to consider receptivity on the upper and lower surfaces of an airfoil at a 
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positive angle of attack where the stagnation point lies on the lower surface, it is more 
convenient to consider w > 0 throughout and take fi > 0 for the upper surface and fi < 0 
for the lower surface. 

4.5.1 Mean flow 

The development of the mean flow is controlled by the mean pressure gradient parameter, 

1 - fiW + fF 

1 + (w - p) 2 

The pressure gradient parameter is plotted as a function of w for three values of p in 
figure 4.2. For p = 0, the pressure gradient is everywhere favorable, falling off mono- 
tonically away from its stagnation point value (/? = 1) and approaching zero as w~ 2 far 
downstream. The behavior for negative w (the lower surface) is identical to that for 
positive w. When aerodynamic loading is introduced, the magnitude of the favorable 
pressure gradient along the upper surface is increased in the region from the stagnation 
point ( w = 0) to the nose (w = p). The pressure gradient parameter then falls rapidly 
to zero at w = p + /z -1 , and the pressure gradient is adverse for all locations farther 
downstream. The adverse pressure gradient decays to zero more slowly with downstream 
distance, becoming proportional to w~ l far downstream. In contrast to the complicated 
behavior of the pressure gradient on the upper surface in the presence of aerodynamic 
loading, the pressure gradient on the lower surface ( w < 0) remains everywhere favorable. 
However, when p is non-zero the pressure gradient on the lower surface decays only as 
M -1 for large distances, and the magnitude of the pressure gradient increases with p. 

As the aerodynamic loading parameter is increased, the position of the minimum wall 
shear moves towards the nose (figure 4.3a), at the same time as the magnitude of the 
minimum wall shear decreases (figure 4.3b). For steady two-dimensional flows, zero wall 
shear corresponds to boundary-layer separation. From Ruban (1982) and Stewartson, 
Smith & Kaups (1982), the minimum wall shear close to the critical angle of attack is 
given by 

4> vn ( w m,V = 0) oc {fl c - fi)*. 

Fitting this behavior to the curve in figure 4.3(b) suggests that separation occurs at 
fi c « 1.15. More accurate determination of separation criteria would require an improved 
numerical scheme, since the current method breaks down near the onset of flow reversal. 

The coefficient Bi(p), which appears in the functions 7)3(77) and p^(rj) arising in the 
large-io asymptotic form of the mean flow (4.29), is undetermined by the asymptotic 
analysis. To determine #1(7/), the wall shear 0 w (io,O) predicted by the asymptotic ex- 
pansion (4.29) is compared with the wall shear obtained by numerical solution of (4.23), 
for each value of fi considered. In fact it proves necessary to modify (4.23) and solve 
for 0(u), 77) — F(tj) in order to obtain the required accuracy when extrapolating for B\ 
in the large - w limit. In figure 4.4, it is seen that B t (0) « 2.08 and that B\{p) increases 
approximately linearly with /j, for —0.8 < fi < 0.8, before decreasing for fi > 0.9. The 
value of B\ is related to the ‘virtual origin’ of the boundary layer, as can be seen by 
noting that the function r/F' — F which B\ multiplies in 7)3 corresponds to the first eigen- 
solution of Libby &; Fox (1963). Their first eigensolution is the streamwise derivative of 
the streamfunction for the Blasius boundary layer. The value i?i(0) relates to a shift of 
virtual origin due to the favorable pressure gradient on the nose of the parabola in the 


(4.46) 
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Figure 4.2: Streamwise variation of the mean pressure gradient parameter, p(w ), for 
H = 0.0, 0.5, 1.0. 
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Figure 4.3: Effect of aerodynamic loading parameter on the wall shear: (a) Position of 
minimum wall shear, w m \ (b) Magnitude of minimum wall shear, 0^(w m ,r; = 0). 
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Figure 4.4: Computed values of the coefficient appearing in the expansion of the 

base flow far downstream (4.29). 

absence of aerodynamic loading; the dependence on fi relates to an additional shift of 
the virtual origin in the presence of aerodynamic loading. 

4.5.2 Unsteady flow 

The value of the Receptivity Coefficient Cy(S, fx), for a particular free-stream disturbance, 
is obtained by comparing numerical solutions of the LUBLE to the asymptotic forms 
(4.34) and (4.35). 

The method used to obtain numerical values for C 3 and C a is described in chapter 2. 
For real w, the real part of t/ 0) (iu) is negative. Thus, the asymptotic eigenfunctions (4.33) 
decay exponentially with downstream distance, while the particular solution ip P remains 
0(1). In addition, the eigenfunctions are inverse ordered, so that the function ^whose 
coefficient is desired decays exponentially faster than all the other eigenfunctions. To 
circumvent these difficulties, the streamwise integration is extended into the complex 
plane in order to make the first asymptotic eigenfunction dominant. Choosing —bn/\2 < 
arg w < — 7r/12, the first eigenfunction becomes exponentially large compared to the other 
components of the unsteady boundary layer solution. However, the pressure gradient 
parameter f3, defined by (4.24), has first-order poles at w = fx±i. Thus for larger values 
of //, arg w must be chosen with care in order to avoid an integration path which passes 
close to the singularity at w = jx — i. For all the values of n considered, integration 
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was performed for arg w = — 0. 1 7r and arg w = — 0.157T, while for smaller values of /x, 
results for arg w = —0.2n and arg w = —0.25/ r were also obtained. For these choices 
of arg w, the solution ip grows exponentially with |iu|. To avoid difficulties associated 
with the rapid growth of ip as a function of |u/|, in our numerical scheme we solve for 
/ = ipexp(-T^) rather than ip. This allows the numerical solution to be continued 
much farther downstream while retaining accuracy (Heinrich Kerschen, 1989). 

Using the wall shear (4.34) as the basis for comparison, we define 


C(/i t S t w) = ^ °) ) (4.47) 

Jo 

where Jq = 0.4356(1 + i ) (5u; 2 /2) ri+1 . Since the argument of w was chosen so that the 
first asymptotic eigenfunction dominates other components of the solution when w » 1, 
/ ~ Ciipi exp (— Ti ) and hence 


C~C: 


1 + Ci 
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W 0 - 774 


+ C2 


In 2 w 
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In w 1 

+ c 3 1- c 4 — I- . . . 

w w 


as w — > oo. (4.48) 


The value of Ci(/x, S) is then obtained by extrapolation. However, for the larger values 
of /x and S considered, a small error in the computed value of B\ leads to a large error in 
C(fi,S,w) for w 1, and hence limits the accuracy of the extrapolated value of C\. In 
the analysis of chapter 2 for a symmetric airfoil in the absence of aerodynamic loading 
(fj, = 0), arg w = — 0.257T was used throughout. With this particular choice, any error in 
the extrapolated value of |Ci| due to inaccuracy in the numerically determined value of 
B\ is eliminated. Clearly this advantage is lost in the present study. However, C\ should 
be independent of the value of arg w used in the numerical integration. By using two 
different values of arg uu in the present analysis, any significant error due to inaccurate 
determination of B i should be readily identifiable. Unfortunately, for the larger values 
of /x and S considered, there is still uncertainty in the numerical value of C\ obtained by 
extrapolation. For fi > 0.9, the relative error may be as large as 25%. 

In figure 4.5, the moduli of the receptivity coefficients, \C 3 \ and |C a |, for free-stream 
disturbances symmetric and anti-symmetric about the nose of the body, respectively, are 
plotted as a function of the aerodynamic loading parameter /x. Results are presented 
for two Strouhal numbers, S = 0.1 and S = 0.3. The overall patterns of the variation 
with the aerodynamic loading parameter /x are qualitatively similar for \C a \ and |C a |. 
However, the receptivity levels for anti-symmetric forcing are somewhat larger than those 
for symmetric forcing. As the aerodynamic loading parameter is increased from zero, the 
receptivity level on the upper surface (/x > 0) decreases until /x « 0.7. The receptivity 
then increases to a local maximum in the vicinity of /x « 0.9, with a subsequent decrease 
in receptivity level as /x is increased further. The quantitative accuracy of the receptivity 
coefficient for values of /x beyond 0.9 is uncertain, however. As explained earlier, results 
close to the critical condition for boundary layer separation are prone to inaccuracies, 
due to difficulties in the extrapolation to large values of w, especially at higher values S. 
A further difficulty is the long downstream distance required for the boundary layer to 
recover from the region of low wall shear induced by the strong adverse pressure gradient 
just downstream of the nose. 

The variations in both |C S | and \C a \ as a function of /x are much larger at S = 0.3 than 
at S = 0.1. For S = 0.3, the receptivity levels on the upper surface drop significantly 
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Figure 4.5: Variation of the magnitude of the receptivity coefficient as a function of 
the aerodynamic loading parameter fi, for two Strouhal numbers, S = 0.1 and 0.3. 
(a) Receptivity due to the free-stream disturbance component symmetric about the 
nose, | Cal; (b) receptivity due to the free-stream disturbance component anti-symmetric 
about the nose, |C 0 |. 


65 



for modest values of aerodynamic loading. Both |C 3 | and \C a \ are less than 0.02 in the 
region 0.2 < /x < 0.6. (It is not possible to provide precise values in this region, because 
the errors that arise in extrapolation are absolute rather than relative.) The subsequent 
rise in receptivity level near /x « 0.9 is particularly striking for the case of anti-symmetric 
forcing with S = 0.3, where \C a \ exceeds 4. 

On the lower surface of the airfoil ( /x < 0), the receptivity coefficients increase as 
aerodynamic loading is introduced, to a global maximum in the vicinity of /x ~ —0.4. 
The receptivity coefficients then decrease toward zero as n approaches —1. For negative 
/x, the value of the receptivity coefficient can be extracted from the computational results 
with a good degree of accuracy, since the pressure gradient is everywhere favorable and 
decays monotonically to zero, with no region of low wall shear. 

The presence of aerodynamic loading adds significant complexity to the development 
of the unsteady flow, especially on the upper surface where the boundary layer first 
experiences a rapid acceleration, followed by a rapid deceleration and then a slow recovery. 
Thus it is difficult to provide a simple explanation for the dependence of the receptivity 
coefficient on the aerodynamic loading parameter. For the symmetric mean flow case 
(/j, = 0) considered in chapter 2, the mean pressure gradient is everywhere favorable, 
monotonically decreasing from its maximum value at the stagnation point to zero far 
downstream as seen in figure 4.2. An increase of the leading-edge nose radius (or Strouhal 
number) was found to decrease the receptivity coefficient. As a speculative explanation 
of this behavior, it was suggested that a favorable pressure gradient in the region where 
the receptivity occurs tends to reduce the receptivity level. At larger values of S, the 
region £ = 0(1) where the receptivity occurs is concentrated closer to the airfoil nose, in 
the region of higher favorable pressure gradient. Hence, the decrease in receptivity for 
increasing S found in chapter 2 is consistent with this explanation. 

However, certain features of the results presented in figure 4.5 are inconsistent with 
this simple explanation. One feature which cannot be explained in this way is the large 
positive peak in receptivity level which occurs on the lower surface at /x w —0.4. It can 
be seen from figure 4.2 that the pressure gradient is negative everywhere on the lower 
surface, and that except for a small region very near the stagnation point, the magnitude 
of the pressure gradient increases monotonically with |/x| . Thus, if the speculative 
explanation discussed above were valid, one would expect the receptivity coefficient to 
decrease monotonically with |/x| in the region /x < 0. A second feature which requires 
further explanation is the rise in receptivity level in the vicinity of /x « 0.9, followed by 
a subsequent decrease for larger values of /x. We address this latter feature first. 

As discussed in §4.1, receptivity comes about through the interaction of free-stream 
disturbances with short-scale streamwise gradients in the boundary layer. In the leading- 
edge mechanism, the receptivity is concentrated in a region of length 0(U e /u> ) extending 
downstream from the mean flow stagnation point. However, for values of \i approaching 
the critical value for the onset of separation, /x c « 1.15, additional short-scale streamwise 
variations develop in the mean flow in the vicinity of the location of minimum wall shear, 
w = w m . From figure 4.3, it is seen that w m « 5 for values of /x > 0.8. The corresponding 
distance from the mean flow stagnation point is s m ~ 15r„, so that us m /U e ~ 15S. Thus, 
for low values of S, the minimum wall shear point lies within the leading-edge receptivity 
region, and the receptivity which occurs in the vicinity of w m can be calculated by the 
present method, while at higher values of S the minimum wall shear point lies downstream 
of the leading-edge receptivity region. 
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The marginal separation receptivity analysis of Goldstein et al(1987) assumes fj. c ~n — 
O(e's') and S = 0(e _1 ). The scaling with respect to /z leads to an interactive viscous- 
inviscid structure for the mean flow in the vicinity of w m , but with a local streamwise scale 
of 0(e 2 ) as compared to the 0(e 3 ) scale of standard triple-deck theory. The fact that the 
length scale of the low wall shear region is larger than that for standard triple-deck theory 
leads to additional complexity relative to the localized receptivity analysis of Goldstein 
(1985), in which the variations in surface geometry and the Tollmien-Schlichting wave- 
length both have the triple-deck scale. Specifically, in Goldstein’s marginal separation 
analysis, the instability waves in the marginal separation region are interactive but their 
wavelength is shorter than the local scale for the mean flow'. This leads to an exponen- 
tially small initial amplitude for the instability wave, but the wave undergoes exponential 
amplification in the marginal separation region and somewhat downstream, so that the 
marginal separation receptivity mechanism may well be important in applications. 

Our analysis assumes a high Reynolds number e«l with /z c — n = 0(1), so that the 
mean flow remains non-interactive in the low wall shear region. Thus, we cannot make 
direct comparisons with the analysis of Goldstein et al (1987). However, as discussed 
above, at sufficiently low Strouhal numbers the region of low wall shear is contained 
within the leading-edge receptivity region, which is treated by our theory. It is then 
natural to investigate the relative importance of the receptivity in the region of low wall 
shear relative to the receptivity near the mean flow stagnation point. 

We investigate the relative importance of receptivity near the mean flow stagnation 
point and receptivity in the region of low wall shear by artificially switching off the 
unsteady slip velocity at various downstream positions, using a smooth transition func- 
tion. Here we are focusing on the relative receptivity for different free-stream disturbance 
functions, and this can be calculated very accurately even when n is large. Essentially, 
the development of the unsteady disturbance along a ray in the complex w plane con- 
tains two phases. First, over a relatively short distance the asymptotic eigenfunction 
grows exponentially relative to the other components of the solution and becomes the 
dominant component. Second, over a much larger distance the asymptotic eigenfunction 
evolves algebraically toward its large-io form (4.33). Both phases of the development of 
the unsteady field must be calculated accurately to determine values of the receptivity 
coefficient. In contrast, a calculation of only the first phase is adequate to determine 
relative receptivity levels, since the same asymptotic behavior of the eigenfunction arises 
in the second phase independent of the form of the free-stream disturbance. 

Thus consider a modified slip velocity 


u 3 (w) = u s (w)$(w,w 0 ) 

where the smooth transition function $(w,w 0 ) is defined by 


<f>(u/,u;o) = 


1 — tanh(z) 

1 — tanh(zo) ’ 


w — w 0 w 0 

z — 7 J Zo — c • 


(4.49) 


(4.50) 


The modified slip velocity is the same as the physical slip velocity at the steady stagnation 
point ( w = 0), then decreases smoothly and monotonically to zero in the vicinity of 
w — w 0 . The transition occurs over a ui-scale of 5 W . The corresponding unsteady pressure 
gradient is then given by 


dp d P . wu s 

d w d w w h 


(4.51) 


67 



and we extract the modified Receptivity Coefficients C 3 (/j,,S,w 0 ,S w ) and C a (ii,S,w 0 ,6 w ) 
from the solutions of 

Hi,,) = Si ^iSk 2 (w -Ii) + ** ~ 2w + - *„«,(«, - M ) | , | 

9 9\ ^ (4.52) 

_ + ^ ~ 

respectively. By varying wq, it is hoped that the relative contribution of different stream- 
wise regions to the total receptivity can be determined. However, the introduction of a 
new artificial stream-wise scale (i.e. S w ) clearly provides another possible source of recep- 
tivity. Thus care must be taken in choosing S w , small enough that the different regions of 
receptivity can be differentiated, but not so small that the receptivity associated with the 
rapid change in unsteady slip velocity dominates other contributions to the receptivity. 
Moreover, <f> is singular when w = w 0 — fyr8 w \. Hence our choice of S w is restricted by 
the condition S w > 2w 0 tan(argic)/7r, which ensures that the integration contour does 
not pass through any such singularity. 

To examine the contributions to the receptivity on the upper surface when the body 
is close to the critical angle-of-attack, in figure 4.6 we plot the magnitude of 


r(n,S,w 0 ,8 w ) 


Cj(A b ^Oi 6 W ) 

cjj^sj 


(4.53) 


for [i = 1.0 and 5 W = 3.0. The results are plotted as a function of the filter switch-off 
point wo, for two Strouhal numbers S = 0.1 and 0.3. The value of |r| is seen to approach 
one as Wq increases, indicating that the receptivity process is complete. The receptivity 
region would be expected to be concentrated closer to the mean flow stagnation point at 
the higher frequency S = 0.3, and this behavior is indeed seen in the plot of |r,|. Upon 
closer inspection, a similar trend can be seen in the plot for |r a |. The minimum wall 
shear point w = 4.54 is marked with a vertical dashed line on the plots. For S = 0.3, the 
receptivity process appears to have been completed by the time the region of low wall 
shear is reached. In contrast, for S = 0.1, |r 3 | does not approach the value of one until 
Wo is approximately 7. Thus, in this case the region of low wall shear appears to make 
some contribution to the total receptivity. This is not unexpected, since the region of 
receptivity us/U e = 0(1) extends farther downstream by a factor of three for the case 
S = 0.1 as compared to the case S = 0.3. However, although the receptivity extends out 
through the region of low wall shear for the case S = 0.1, there is no indication in figure 
4.6(a) of a concentration of receptivity in the region of low wall shear. It is interesting 
to note that |r a | is not a monotonically increasing function of wq. This feature indicates 
that, at least for |r s |, large receptivity close to the leading edge is partially cancelled 
out by receptivity processes farther downstream. We should note that some caution is 
necessary in evaluating the results for S = 0.1, since in that case the wavelength of the 
unsteady motion is longer and the results could have been influenced more by the finite 
value of the filter width 6 W . 

The results presented in figure 4.6 show that the receptivity region extends farther 
dowmstream at lower values of the Strouhal number. In order to examine the extent 
to which it is possible to scale this feature out of the results, in figure 4.7 we plot |F a | 
as a function of S^Wq, for a constant filter width in terms of the receptivity variable, 
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Figure 4.6: Plot of |T| for ^ = 1.0, 8 W = 3.0, and S = 0.1 (solid line) and S = 0.3 (dashed 
line), illustrating the relative receptivity as a function of the streamwise location of the 
filter, w 0 , for (a) component of the disturbance symmetric about the nose; (b) component 
of the disturbance anti-symmetric about the nose. The dotted line marks the position of 
minimum wall shear, w = 4.54. 


= S$6 W = 0.96. Results for presented for two values of the Strouhal number, 5 = 0.1 
and 0.3. The case /x = 0 is shown in figure 4.7(a). Here, plotting the results in terms of 
a variable scaled by the receptivity length scale has produced an excellent collapse of the 
results for the two values of S. In figure 4.7(b) for the case /x = 0.5, the collapse is not as 
good, with the curve for S = 0.3 extending out to a somewhat larger value of before 
approaching the value of one. This trend which appeared in figure 4.7(b) is magnified 
further in figure 4.7(c) where results for /x = 1.0 are presented. However, in all cases 
the leading-edge receptivity process is completed by the time a downstream distance of 
S*w 0 = 3 is reached. A number of factors may contribute to the lack of collapse in figures 
4.7(b,c). First, the flow along the upper surface of the airfoil becomes successively more 
complex as the aerodynamic loading parameter is raised, so that there is not even an 
approximate similarity behavior for the mean boundary layer whose streamwise gradients 
are responsible for the receptivity process. Second, the strong acceleration of the flow 
between the stagnation point and the nose of the airfoil may introduce a significant shift in 
the ‘virtual origin’ of the boundary layer which could affect the streamwise development 
of the receptivity process. Thus, it seems unlikely that a collapse of the results for 
different frequencies could be achieved for cases with significant aerodynamic loading. 

Similar calculations were performed for /x = —0.4 and —0.8. The results for |T| 
indicate that the receptivity on the lower surface is also concentrated in the region S 2 ic 0 < 
3. The calculations did not produce any features which provide an explanation for the 
high receptivity levels in the vicinity of /x = —0.4. 

Despite the fact that not all aspects of the results are well understood, a number of 
general conclusions can be drawn. It is clear that modest levels of aerodynamic loading in 
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the leading-edge region cause a decrease in the receptivity level for the boundary layer on 
the upper surface of the airfoil, and an increase in receptivity for the lower surface. The 
effects are more pronounced at higher values of the Strouhal number, where the region 
of receptivity is concentrated nearer the stagnation point. Thus, it seems reasonable to 
suppose that the effects are related to the behavior of the mean flow pressure gradient 
near the stagnation point. For the upper surface, the introduction of aerodynamic loading 
leads to a movement of the stagnation point toward the lower surface, and an increased 
favorable pressure gradient in the region between the stagnation point ( w = 0) and the 
leading edge ( w = //). The decrease in receptivity level associated with these effects is 
consistent with the results obtained in chapter 2 for the case fx = 0. The decrease in 
receptivity level on the upper surface in the presence of modest aerodynamic loading 
is of significance for applications such as laminar flow design. The subsequent rise in 
receptivity for values of fx nearer the critical value for separation also has important 
practical implications. For the lower surface, an argument that the increase in receptivity 
with modest aerodynamic loading is related to the pressure gradient near the stagnation 
point appears much more tenuous, since the changes in the pressure gradient near the 
stagnation point are not pronounced. However, from a practical standpoint, receptivity 
levels for the lower surface are usually of less interest, since the pressure gradient on the 
lower surface is typically favorable and the instability waves have much smaller growth 
rates. 

It must also be kept in mind that, for leading-edge receptivity, quantities such as the 
instability wave amplitude at the neutral stability point depend crucially on a combi- 
nation of receptivity and stability properties. There can be a subtle interplay between 
these, since a small change in growth/decay rate can counteract a significant change in 
receptivity level, for example. In fact, for the leading-edge receptivity mechanism, even 
the issue of the most appropriate definition of receptivity level is not straightforward. 
Experimental and computational researchers often favor a definition based on the lower 
branch amplitude. This is certainly of practical utility and has advantages with respect 
to measurements, but it does not distinguish between the receptivity and stability as- 
pects of the disturbance development. The definition we have adopted has the advantage 
of separating the receptivity aspects from stability considerations, but of course the re- 
sults must be combined with a stability calculation in order to predict amplitudes in the 
Orr-Sommerfeld region farther downstream. 
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Chapter 5 
Conclusions 


In this report, we have examined the influences of airfoil thickness and mean aerodynamic 
loading on leading-edge receptivity to acoustic free-stream disturbances. A high Reynolds 
number asymptotic analysis (e = ( uu/U 2 )~ e < 1) has been developed for the case of a 
two-dimensional airfoil with a parabolic leading edge. The airfoil is assumed thin (6 1) 

and the mean flow Mach number M = U/c is assumed small. The receptivity process 
takes place in the region of the boundary layer where the distance from the stagnation 
point is 0(U/uj). Previous analyses of leading-edge receptivity considered the case of the 
Blasius boundary layer on a semi-infinite flat plate. In the present analysis, we assume 
that the radius of curvature r„ of the airfoil nose is 0(U/uj). Thus, the leading-edge 
receptivity process takes place in the vicinity of the airfoil nose, and the receptivity 
coefficient depends on the thickness distribution of the airfoil only through the Strouhal 
number, S = u>r n /U , which has 0(1) values. 

For a symmetric airfoil at zero angle of attack, the mean flow is symmetric about 
the leading edge and the receptivity characteristics on the two sides of the airfoil are 
the same. In contrast, for an airfoil with camber or non-zero angle of attack, there is 
generally a mean flow component around the leading edge and the boundary layers on 
the two surfaces of the airfoil have quite different characteristics. For S = 0(1), the 
leading-edge receptivity process depends on the mean aerodynamic loading only through 
a parameter related to the local flow in the vicinity of the airfoil nose. The aerodynamic 
loading parameter which enters the leading-edge receptivity analysis is 

H = 2a e {6 2 b/r n )?, (5.1) 

where b is the airfoil semi-chord and a e is an effective angle of attack parameter for 
the leading-edge region, discussed below. An aerodynamically loaded body typically has 
a mean flow component around the leading edge, from the lower surface to the upper, 
corresponding to positive values of a e and //. The parameter fi is a measure of the 
asymmetry of the mean flow in the leading-edge region. Separation of the boundary 
layer, on the upper surface just downstream of the leading edge, occurs when fj. reaches 
the critical value // c ss 1.15. When discussing results below, positive values of /r are 
taken to correspond to receptivity on the upper surface for a typical aerodynamically 
loaded body, while negative values of fi correspond to receptivity on the lower surface for 
a typical aerodynamically loaded body. 

The parameter /i is seen to depend on the airfoil thickness ratio 6, semi-chord b and 
nose radius r„, in addition to the effective angle of attack parameter a e . A brief discussion 
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of these dependencies seems worthwhile. The strength of the mean flow component 
around the leading edge depends on both the geometric angle of attack and the camber 
distribution of the airfoil; these dependencies are contained in the parameter a e , defined 
by (4.8). For a thin airfoil, the camber and angle of attack must be small in order to 
avoid flow separation; hence, in defining a e the small parameter S has been extracted 
from the physical angle of attack and from the camber function, so that a e = 0(1). For 
a thin airfoil with a parabolic nose, the nose radius of curvature r n is 0{5 2 b), so that the 
quantity under the radical is also 0(1). Physically, when the airfoil thickness is decreased 
the nose radius also decreases, leading to higher flow velocities around the sharper edge, 
so that the camber and/or angle of attack must be decreased correspondingly in order 
to avoid flow separation. 

The asymptotic structure for high Reynolds number (e 1) involves two stream- 
wise regions. In the first region, corresponding to downstream distances of 0(U/u), 
the unsteady component of the motion in the boundary layer is governed by the lin- 
earized unsteady boundary layer equation (LUBLE). Here the inviscid pressure field and 
slip velocity induced by the free-stream disturbances drives the unsteady motion in the 
boundary layer. The leading-edge receptivity process takes place in the LUBLE region. 
The second region is farther downstream, at distances of 0(e~ 2 U/u>). In this region the 
unsteady disturbance is governed by the triple-deck structure, corresponding to the high 
Reynolds number asymptotic form of the Orr-Sommerfeld equation (OSE) in the vicin- 
ity of the lower branch. The free-stream disturbance that is present in the OSE region 
generates only a Stokes wave in the boundary layer, which is unrelated to the TS wave 
eigensolution. Thus, the amplitude of the TS wave cannot be determined by analysis of 
the OSE region alone. Rather, the amplitude of the instability wave is determined by 
the asymptotic matching of the LUBLE and OSE regions. 

The asymptotic matching of the LUBLE and OSE regions is carried out in the over- 
lap region, 1 <C ujx * jU -C c -2 , where x* is the dimensional distance from the stagnation 
point. (The analysis was carried out in a parabolic coordinate system, but the physi- 
cal variable x* is utilized in this discussion for convenience.) Thus, the overlap region 
corresponds to distances ‘far downstream’ in the LUBLE region, but ‘far upstream’ in 
the OSE region. Far downstream in the LUBLE region, the unsteady flow consists of 
a generalized Stokes wave driven by the local unsteady pressure gradient in the free 
stream, plus an infinite set of asymptotic eigensolutions, ipi. These asymptotic eigen- 
solutions were first found by Lam &: Rott (1960) for the case of the Blasius boundary 
layer. The asymptotic eigensolutions of the LUBLE obtained in the present work are 
generalizations of the Lam-Rott expression, taking account of the influences of airfoil 
thickness and aerodynamic loading. The Lam-Rott eigensolutions are ‘asymptotic’ in 
the sense that they exist only for ux* jU 1, where the mean flow varies slowly in the 
streamwise direction relative to the scale U/oj. The Lam-Rott eigensolutions contain ve- 
locity but not pressure fluctuations, and are uniformly valid solutions of the LUBLE for 
ux*/U ;» 1. (However, they are not uniformly valid solutions to the full Navier-Stokes 
equations, as discussed below.) The asymptotic analysis for uix*/U » 1 determines the 
form of the Lam-Rott eigensolutions, but not their coefficients C*. These coefficients can 
be obtained only through a full solution of the LUBLE, which must be determined by 
numerical methods. 

The wavelengths of the Lam-Rott eigensolutions shorten progressively with distance 
downstream. When the OSE region is reached ( ux*/U — 0(e -2 )), the self-induced 
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pressure field associated with the displacement thickness of each asymptotic eigensolution 
becomes significant, and the triple-deck structure replaces the LUBLE as the correct 
asymptotic approximation to the Navier-Stokes equation. The asymptotic matching of 
the LUBLE and OSE regions in the overlap domain 1 •C u>x* jU -C e~ 2 shows that the 
first Lam-Rott eigensolution ipi of the LUBLE matches onto the TS wave eigensolution 
of the OSE. Thus, the form of the free-stream disturbance and the geometry close to 
the nose influence the amplitude of the TS wave only through the coefficient C\ of the 
first asymptotic eigensolution. Therefore, we call C\ the ‘receptivity coefficient’. The 
primary objective of this work is to determine the receptivity coefficient as a function 
of the leading-edge nose radius parameter S, aerodynamic loading parameter /i, and 
free-stream disturbance characteristics. 

In order to determine the receptivity coefficient, a numerical solution of the LUBLE is 
obtained for specified values of the nose radius parameter S and the aerodynamic loading 
parameter /z. The free-stream disturbance must also be specified; it can be represented 
as a superposition of components symmetric and anti-symmetric about the nose, as dis- 
cussed below. The value of the receptivity coefficient C\ is then found by comparing 
the numerical solution for large downstream distances ( vx*/U >• 1) with the analytical 
expression for the first asymptotic eigensolution, tpi. However, this procedure cannot be 
carried out for real values of x*, because the first asymptotic eigensolution becomes ex- 
ponentially small far downstream. Moreover, the asymptotic eigensolutions are inversely 
ordered, so that the higher-order eigensolutions ipi, i — 2,3,... are exponentially large 
compared to Tp\. Thus, it is extremely difficult to extract the coefficient C\ from the 
numerical solution for real x*. This difficulty is circumvented by taking advantage of 
the fact that the LUBLE is analytic in x*, and extending the computation into a region 
of the complex x*-plane where the asymptotic eigensolutions grow exponentially with 
downstream distance. The inverse ordering of the eigensolutions is also rectified by this 
process, so that ipi is exponentially dominant far downstream, allowing its coefficient to 
be extracted from the numerical results. A number of analytical features are exploited 
in implementing the numerical computations, in order to increase the precision of the 
results for C\. However, calculations for conditions very near boundary layer separation 
are difficult, so that despite these efforts there is still uncertainty in the numerical value 
of Ci for values of ji near /z c rs 1.15. 

For sufficiently small amplitudes of the unsteady motion, the receptivity and the 
initial growth of the instability wave are linear processes. Our theory addresses this linear 
regime. Thus, it is natural to define the receptivity level as a ratio of the amplitude of 
the instability wave in the boundary layer, or its upstream precursor, to the unsteady 
motion in the free stream. For localized receptivity which occurs in the OSE region, 
an instability wave is generated directly at the receptivity site (e.g. a short-scale wall 
hump), so that the receptivity level can be expressed in terms of the instability wave 
amplitude at that point. In contrast, for leading-edge receptivity, the receptivity process 
occurs upstream in the LUBLE region. In this case the ‘best’ definition of receptivity 
level is not so obvious. Various measures of receptivity can be considered. We have 
chosen to define Cj, the coefficient of the first Lam-Rott asymptotic eigensolution, as 
the ‘receptivity coefficient’. 

Alternative definitions of receptivity level that have been considered in computations 
and experiments are values based on the TS wave amplitude at the lower neutral-stability 
point (branch I), or an extrapolation of the TS wave amplitude back to the leading edge. 
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A receptivity level defined in terms of the TS wave amplitude at the lower branch neutral- 
stability point does allow easy comparison with experiments or computations, but has the 
disadvantage that the results depend on both frequency and Reynolds number, making 
presentation of the results much less compact. Extrapolation of the TS wave amplitude 
back to the leading edge is not advisable, since the TS wave is not a valid solution of 
the disturbance equations near the leading edge. In fact, in extrapolating back to the 
leading edge, the slowly varying amplitude A{x*) multiplying the mode-shape function 
of the TS wave should be considered, and A{x*) is singular as i* -* 0, becoming infinite 
at the leading edge. 

From a theoretical standpoint, the coefficient C\ of the first Lam-Rott asymptotic 
eigensolution of the LUBLE is the most attractive measure of receptivity. This Lam- 
Rott eigensolution evolves into the unstable TS wave farther downstream so that, when 
appropriate scaling factors related to the asymptotic matching of the LUBLE and OSE 
regions are introduced, C\ is also the coefficient of the TS wave. An extremely attrac- 
tive feature of this receptivity measure is that C\ is independent of the frequency and 
the (asymptotically large) Reynolds number, thus providing the receptivity results in 
the simplest dimensionless form. The quantity Cj is determined by the unsteady flow 
behavior upstream of the region of instability and therefore focuses on the receptivity 
aspect of the unsteady flow development, with much less influence from the global stabil- 
ity properties of the flow. Also, for fixed values of the frequency and Reynolds number, 
variations in the receptivity coefficient C\ with changes in the free-stream disturbance 
characteristics translate directly into variations in the TS wave amplitude. 

The instability wave amplitude which is found for a given free-stream disturbance in 
a specific situation depends on a number of processes. First, the free-stream disturbance 
interacts with the body (and other surrounding surfaces, if present) to produce an inviscid 
unsteady pressure gradient and slip velocity which drives the unsteady motion in the 
boundary layer. Next, through the receptivity process, some of the energy contained in 
the unsteady motion just outside the boundary layer is transferred to the first Lam-Rott 
asymptotic eigensolution, which eventually evolves into the TS wave of the OSE region. 
Obviously, both the inviscid processing and the subsequent receptivity process influence 
the amplitude of the TS wave. 

The present work has addressed receptivity to free-stream acoustic waves. The dis- 
turbance velocity of the acoustic wave far from the airfoil has been utilized in normalizing 
the receptivity coefficient C\. The acoustic wave is assumed to propagate in a direction 
perpendicular to the airfoil span, so that the unsteady interaction is two-dimensional. 
The angle between the direction of propagation and the airfoil chordline is denoted by 
9. The Mach number has been assumed small, so that the acoustic wavelength is long 
compared to the hydrodynamic length scale U/u> for the unsteady motion. Thus, in 
analyzing the LUBLE region, the flow can be considered incompressible. 

Since S has been assumed 0(1), only the characteristics of the free-stream disturbance 
on the scale of the nose radius, r„, are required. The incompressible irrotational flow in 
the vicinity of the airfoil nose can be represented by a superposition of two components: 
(a) a symmetric flow past the parabolic nose, and (b) an anti-symmetric flow around 
the parabolic nose. The disturbance slip velocity which drives the unsteady motion 
in the boundary layer can similarly be separated into symmetric and anti-symmetric 
components, with coefficients k s (6) and K a {0), respectively (2.50). Since the equations 
governing the unsteady motion in the LUBLE region are linear, the response to the 
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symmetric and anti-symmetric components of the disturbance slip velocity near the airfoil 
nose can be superposed, so that the receptivity coefficient can be expressed in the form 
(2.52) 

Ci = K„(8)C s (S,fi) + K a (0)C a (S,n). (5.2) 

This expression effectively separates the effects of the inviscid processing (contained in 
k, 3 and k q ) from the intrinsic receptivity of the boundary layer. The quantities C 3 and 
C a are the receptivity coefficients to the symmetric and anti-symmetric components of 
the disturbance slip velocity, respectively. The receptivity coefficients C 3 and C a are 
universal functions which do not depend on the propagation angle of the acoustic wave 
or the mean flow Mach number. Results for C 3 and C a are discussed below. 

In a given physical situation, the behavior observed for the receptivity as a function 
of various parameters may be influenced strongly by the inviscid processing of the dis- 
turbance as represented by the coefficients k 3 and n a . In addition to the incidence angle 
of the acoustic wave, k 3 and n a depend on the acoustic reduced frequency k = ub/c. 
Relatively simple expressions for k 3 and K a can be obtained for the two limiting cases of 
k 1 and k 1, as discussed in section 2.4. The acoustic reduced frequency can be ex- 
pressed as the product k = aM , where a = ub/U is the aerodynamic reduced frequency. 
For extremely low Mach numbers, the acoustic wavelength is long not only compared to 
the hydrodynamic length scale U/u), but also compared to the airfoil chord. Then k 1 
and we obtain (2.58), 

k 3 = cos 6, K a = a* sin0. (5.3) 

Note that a = Sb/r n , which is a large quantity for a thin airfoil with S = 0(1). The 
second limit is that of M small but a very large. Then k 1 and the acoustic wavelength 
is short compared to the airfoil chord (but still long compared to the airfoil nose radius). 
For k ;» 1, we obtain (2.62) 


k 3 = cos 6 , 
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In both limits, the coefficient for the anti-symmetric component of the disturbance is 
large compared to the coefficient for the symmetric component. Basically, the unsteady 
disturbance level near the airfoil nose is much larger for an oblique acoustic wave than 
for a wave propagating parallel to the airfoil chord (6 = 0), due to the induced flow 
around the relatively sharp leading edge of the thin airfoil. For k <§; 1, the size of 
n a is determined by finite chord effects, while for /c 2> 1 the size of K a is determined 
by compressibility effects. It should be noted that the dependence of K a on 9 is quite 
different for the two limits 1: « 1 and k 1. These results are for the case of an 
obliquely propagating acoustic wave in free space. When the acoustic field is contained 
in an enclosed region, such as within a wind tunnel, more complicated behavior is found 
(see Heinrich &; Kerschen 1989). 

In chapter 2, results for C s (S ) and C a (S) are presented for the case of a symmetric 
mean flow past a parabolic leading edge (// = 0). The magnitude of C s is plotted as a 
function of S in figure 2.4. The magnitude of C s at first rises slightly for very small S 
(0 < S < 0.02), then decays monotonically with further increases in S. For S = 0.3, 
the value of |C,| is only about 15% of its value for S = 0. The decay of \C 3 \ appears 
to be exponential in S, although this has not been verified analytically. The phase, 
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arg(C 3 ), also plotted in figure 2.4, increases monotonically with 5. The increase in 
arg(C a ) corresponds to a lag in phase of the TS wave. 

The magnitude and argument of the ratio of receptivity coefficients, A = C a /C 3 , is 
plotted in figure 2.5. The magnitude |A| is approximately 6 in the flat-plate limit 5 — * 0, 
where the anti-symmetric component of the slip velocity has a square-root singularity at 
the leading edge. The value of | A j decreases rapidly as a finite nose radius is introduced, 
relieving the singularity at the leading edge. Thus, the behavior of C a for very small 5 is 
quite different from that for C 3 , where a small rise in magnitude was seen. The behavior 
of C a for 5 -C 1 is in fact singular, as discussed below. Beyond 5 = 0.05 the decline in 
|A| becomes more gradual, the value of |A| remaining above 2.5 right out to 5 = 0.3. For 
all values of 5 examined, the anti-symmetric component of the slip velocity was found 
to be more effective than the symmetric component in generating a TS wave. 

The important role of the inviscid processing of the free-stream disturbance is illus- 
trated in figures 2.6 and 2.7. In figure 2.6, the receptivity coefficient C\ is plotted as a 
function of the incidence angle 6 of the acoustic wave, for the case 1 and a = 10. 
The overall pattern is similar to the function |sin0|, indicating that the dominant con- 
tribution comes from the receptivity to the anti-symmetric part of the disturbance. This 
behavior is explained by noting that the ratios |/c a /K; a | and |A| = \C a /C 3 \ are both sig- 
nificantly greater than one. In figure 2.7, the corresponding results are plotted for the 
case k 1 and M = 0.1. The overall pattern in this case is similar to the function 
| sin j0|, indicating that the dominant contribution again comes from the receptivity to 
the anti-symmetric part of the disturbance, since |« a //c 3 | is significantly greater than one 
in this limit as well. 

The results presented in figure 2.5 suggest that C a is singular in the limit 5 — ► 0. In 
chapter 3, a small-5 asymptotic theory is developed for the case of a symmetric mean 
flow past a parabolic leading edge (/i = 0). For small values of 5, the nose radius 
r n of the airfoil is small compared to the streamwise length scale U/u for the LUBLE 
region where the receptivity occurs. Thus, in this case three streamwise regions must 
be considered: the nose region (which scales on r n ), the LUBLE region (which scales on 
U/co) and the OSE region (which scales on e~ 2 U/u>). The solution in the nose region is 
a quasi-steady perturbation at leading order. In contrast, the solution in the LUBLE 
region is fully unsteady, and viscous-inviscid interactions appear in the OSE region. 
The matching of the LUBLE and OSE regions follows that discussed above. From the 
forms of the expansions in the nose and LUBLE regions, the receptivity coefficient to 
the symmetric component of the free-stream disturbance, C 3 (S), has an expansion in 
terms of 5°, 5 log 5, 5, 5 1887 , ...; the corresponding expansion for C a (S ) contains 5°, 
55, 5 log 5, 5, 5 1 ' 387 , .... However, asymptotic matching of the (numerical) solutions 
for the nose and LUBLE regions leads to the surprising result that the coefficients of the 
5 log 5 terms in C 3 and C a are approximately zero (i.e. to within the accuracy of the 
computations). The following asymptotic expansions are obtained: 

C,(5) ~ (-0.441 + 0.841 i) + (0.62 + 3.67 i)5 + 0(5 1887 ), 

C a (S) ~ (-5.33 + 1.66 i) 4- (12.08 + 6.34 i)5* + (-10.7 + 24.2 i)5 + 0(5 1387 ). 

These small-5 asymptotic results are compared to numerical results in chapter 3. The 
small-5 expansions are found to have only a small range of validity near 5 = 0. However, 
the singular structure of C a (S) has been elucidated, and the asymptotic formula provides 
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accurate quantitative results in the region near 5 — 0 where the receptivity coefficient 
varies most rapidly with respect to this parameter. 

The influence of aerodynamic loading on the leading-edge receptivity process was 
investigated in chapter 4. The analysis of the Lam-Rott asymptotic eigensolution is 
much more complicated in this case, due to the presence of two parameters (5 and /j) 
and the need to retain many more terms in the asymptotic expansions. The numerical 
computations are also more demanding. For = 0, a single computation is adequate to 
fully describe the mean flow in the vicinity of the parabolic leading edge. In contrast, in 
the presence of aerodynamic loading, a separate computation for the mean flow is required 
for each value of n. For conditions in the vicinity of boundary layer separation ( fi c ~ 1.15), 
the mean flow calculation becomes much more difficult since the behavior of the boundary 
layer becomes very sensitive to small changes in parameter values, and a large downstream 
distance is required for the boundary layer to recover from the region of low wall shear. 
The numerical extraction of the receptivity coefficient becomes correspondingly more 
difficult. Thus, for conditions near boundary layer separation, values of the receptivity 
coefficient cannot be obtained to the same accuracy as elsewhere. 

In figure 4.5, results for |C S | and \C a \ are presented as a function of the aerodynamic 
loading parameter /x, for two values of the Strouhal number, 5 = 0.1 and 0.3. The overall 
patterns of the variation with n are qualitatively similar for |C,| and \C a \- However, the 
receptivity levels for anti-symmetric forcing are somewhat larger than those for symmetric 
forcing, as found previously in the absence of aerodynamic loading (see figure 2.5). As 
the aerodynamic loading parameter is increased from zero, the receptivity level on the 
upper surface (fi > 0) decreases until /i « 0.7. The receptivity level then increases to a 
local maximum in the vicinity of n « 0.9, with a subsequent decrease in receptivity level 
as fj. is increased further. On the lower surface (/i < 0), the receptivity level increases 
as aerodynamic loading is introduced, to a global maximum in the vicinity of n « —0.4. 
The receptivity level then decreases as [i approaches —1. The variations in both \C S \ and 
\C a \ as a function of n axe much larger for 5 = 0.3 than for 5 = 0.1. 

As discussed earlier, receptivity can occur in the leading-edge region, or in a localized 
region farther downstream in the boundary layer where some external feature forces the 
boundary layer to adjust on a short streamwise length scale. For conditions near bound- 
ary layer separation where fx c — (j. = O(e^), the mean boundary layer develops a short 
local streamwise length scale of 0(e 2 ). Goldstein, Leib & Cowley (1987) developed a local 
receptivity analysis for this ‘marginal separation’ case, also assuming that 5 = 0(e -1 ). 
The streamwise length scale of 0(e 2 ) for marginal separation is longer than the standard 
triple-deck length scale of 0(e 3 ) associated with localized receptivity, leading to addi- 
tional complexity in the analysis. We cannot make direct comparisons with the results 
of Goldstein et al., since our analysis assumes — n = 0(1) and 5 = 0(1). However, 
at sufficiently low Strouhal numbers, the region of low wall shear is contained within 
the leading-edge receptivity region and hence is treated by our theory. We investigated 
the relative importance of receptivity in the region of low wall shear for these cases by 
artificially switching off the unsteady slip velocity at various downstream positions, using 
a smooth transition function. For small values of 5 and larger values of /i, within the pa- 
rameter range examined, the region of low wall shear appears to make some contribution 
to the receptivity, but it is not dominant over the receptivity from the region nearer the 
leading edge. It should be noted that this observation does not contradict the conclusions 
of Goldstein et al., since they addressed a different region of parameter space. 
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In summary, the leading-edge receptivity theory has been extended to address the 
influences of leading-edge thickness and mean aerodynamic loading. Receptivity coef- 
ficients have been calculated for the case of acoustic free-stream disturbances in a low 
Mach number flow. A number of general conclusions can be drawn. The introduction of 
leading-edge thickness causes a decrease in receptivity levels relative to the case of a fiat 
plate. It is also clear that modest levels of aerodynamic loading decrease the receptivity 
on the upper surface, but increase the receptivity on the lower surface. The effects are 
more pronounced at higher values of the Strouhal number, where the region of receptiv- 
ity is concentrated nearer the stagnation point. A rise in receptivity level on the upper 
surface occurs for larger values of aerodynamic loading. The decrease in receptivity level 
on the upper surface in the presence of modest aerodynamic loading is of significance for 
applications such as laminar flow design. The subsequent rise in receptivity for values of 
/i nearer the critical value for separation also has important practical implications. 

It must be kept in mind that, for leading-edge receptivity, quantities such as the 
instability wave amplitude at the lower branch neutral-stability point depend crucially 
on a combination of receptivity and stability properties. There can be a subtle interplay 
between these, since a small change in growth/decay rate can counteract a significant 
change in receptivity level, for example. There is a need to relate receptivity levels defined 
in terms of C\ to levels expressed in terms of branch I amplitudes. In the asymptotic 
theory, this is achieved by developing the asymptotic expression for the TS wave in 
the OSE region and matching this expression to the generalized Lam-Rott asymptotic 
eigensolution. An alternative approach is the numerical solution of the disturbances 
equations, using the Lam-Rott asymptotic eigensolution as a starting condition and 
continuing downstream of branch I. 

Important questions for further study include the receptivity to vortical free-stream 
disturbances, and compressibility effects for high subsonic and supersonic Mach numbers. 
Corroboration of the present theory with experimental and computational results is also 
an important objective. In comparing results of the asymptotic theory with computations 
or experiments, it is essential to take proper account of differences in the definitions of 
receptivity levels. 
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Appendix A 

Evaluation of integrals involving 
Airy functions 


Here we obtain alternative expressions for the integrals Ji, defined by 

TOO 

Ji= [z - 2 o]*Ai 2 (z) dz, (A.l) 

J zq 

where Ai'(z 0 ) = 0. Setting w — Ai 2 (z), this satisfies 

vJ" — 4 zw' — 2 w = 0 (A. 2) 

(Abramowitz & Stegun (1965); equation 10.4.57), and hence 

J w(z') d z' = — — Azw). (A.3) 

Thus Jq{zq) follows, noting that ui"(zq) = 2zoAi 2 (zo). The other integrals Ji are then 
obtained by successive integration by parts, 

1 r°° 

Ji = - / (z - z 0 y~ l (w" - Azw ) dz 

2 JZQ 

1 r°° 

= —2iJi — 2izoJi_i + - (z — Zq)'~ 1 w" dz. (A. 4) 

2 JzQ 


Hence, 


Jo 

Ji 

Ji 

J h 


= -z 0 A\ 2 (zo), 

3 Z 0 Jo j 

= g( AzqJi + Ai 2 (zo)), 

2i _ T *'(*'- !)(* - 2 ) j 

2 oA-i H „/, . „-s — Ji- 3, 


1 +2 i 


2(1 +2i) 


i > 3. 


(A.5) 
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Appendix B 

Matching of asymptotic 
eigenfunctions to 
Tollmien-Schlichting modes 


As the wavelength of the asymptotic eigenfunctions of the Linearized Unsteady Boundary 
Layer Equations progressively shortens with distance downstream, an outer inviscid layer 
begins to interact with the wall viscous layer. This occurs when £ = 0(e _l ). Writing 
£ = ef, we anticipate the rapid variation of ip' by writing 

i>' = G{r],£)exp(ijkd£/e 3 ) (B.l) 

where the scaling is dictated by the need to match back to (3.17) as £ — ♦ 0. As in the 
LUBLE regime, the exponential streamwise variation /c(£) is determined as an eigenvalue 
from the matching between the different layers. In the Orr-Sommerfeld region, the main 
deck (77 = 0 ( 1 )) is matched to an outer irrotational layer (77 = 0 (e -1 )) and to the inner 
viscous layer (77 = 0(e)). When £ = 0(e -1 ), the classic triple deck scalings arise and 
the effect of the outer layer enters the expression for k at leading order. In the present 
investigation, fc(f) was determined using direct matching of the stream function between 
the decks, rather than by trying to modify standard Orr-Sommerfeld results to include 
wall curvature and mean pressure gradients effects. 

Details of the matching are not included here, but an implicit expression for k(£, e) 
is eventually obtained, 


-^e 2 ) 

e ) 


Q ( z 0) = 7 3 - h ( z 0 , e , i ), 


(B.2) 


where 


f)t 7 \ - e^oAi'fyo) _ ( —iff 6 V 

Q{ o) - /» Ai(S)d 5 ’ 20 UVI 2 * 2 ) ’ 


7 = 1 + A X ^ SJP - + B - 1 


£ 2 /S ?/S 


and 


h{z 0 ,e,i) = | 3 + e i,r/4 | 2 Zo /2 (? - e + |z^i fl + e 2 + 0(e 3 lne), 

\ 12 0 / \ 1Z 0 z o J 
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Here £ = £/£/q and J* are numerical constants involving integrals of the mean flow through 
the boundary layer, given by Goldstein (1983: equations 4.44-4.45). Solving (B.2) for 
zo(Cc) 6 ives 


Zo — Zq 



fr \ 

zoQ'(zo)/ 


(B.3) 


where z 0 ~ zoo + ezoi 4- e 2 zo2 + 0(e 3 In e) is the flat-plate result, and 


T = —3A 


?/S ^ B + 2 \ 2 /S °( e2lnc )- 


(B.4) 


Here we have written 7 and T in terms of £ rather than £ for conciseness. Thus we finally 
obtain 


k = k\l-A 


in e/s 
e/s 


+ (§ ~B) 


e/s, 


1 - 


3£ 3 r 


‘ZZqQ^Zq ), 


(B.5) 


where k = e 3,7r//4 ^ 2 2 0 3/,2 /C/q is the flat plate result. 

The leading order term in the small-e expansion of z 0 is given implicitly by 


Q{zq 0) = | 3 . 


(B.6) 


which has a set of roots Zoo(|). As £ — > 0, Zqq — > —pi, where p { is defined in §3, namely 
the i-th root of Ai'(— p) = 0. The set of roots k^(£) then follows from (B.5). As £ — > 0, 


Jfe® 


~ e 


, e *'*£ 2 


3/2 


\- A ±?lS + (1 
v A £ 2 /s +(a 


B) 


?/S, 


and hence 




e _ 4 i,r ^ 3 

Uifif 


5 - + (2 A, + 3/2 - B,)^ 


S' 

■f 2 


(B.7) 


(B 8) 


which matches to the exponential variation of the eigenfunctions of the LUBLE region 
(3.17). Numerical solution of Q(zq 0 ) = £ 3 shows that as £ increases, the imaginary part 
of the leading order term in fcW stays positive for i = 2,3, . . ., corresponding to damped 
modes, but that lm(k^) < 0 for £ > 3.03, which corresponds to a growing wave. Thus 
we have demonstrated that the first LUBLE eigenfunction does indeed match on to the 
unstable Tollmien-Schlichting mode. Moreover, (B.5) can be used to calculate the change 
in the position of the neutral stability point compared to the flat plate case. 
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Appendix C 

/ • \ 

Definition of the operators Afp iy(w) 
and TZp^w) 


In §§3.3 and 3.4 a set of partial differential operators Afp^w) are used. These are defined 
as 


O) = 
O) = 

Af^{w) = 


w Tm + F w w + [*£ 2 + (0 4- \£ 3 )F']wjj + (7 — A £ 3 )F"w ^ 
+Z(F"w< - F^t), 


where 


(C.l) 


#{?» = G^^(/3 + AC 3 )GX-(^-AaG>-C(G> ? -GX ? ), ' 
~ (2Gi — G 2 )uv, + ((/3 + Af 3 )G£ — 2(Fjf- Gj ))«;,, 

+(7 ~ A^ 3 )G> 4- f (GJ^ — G^Wtf) 


(C.2) 


In addition, the set of differential operators Ttp\(f) are used, defined by 


*£(/) = Gi/„ - 7Gj7 1 , c 3 . 

Rgi(/) = (2G. - G 2 )/„ + (|Sff 2 - 2(F + G',))/, + idf j 

so that 

(C.4) 

for i = 3, 4. 
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Appendix D 

Asymptotic form of base flow far 
downstream 


In order to determine the form of the unsteady perturbation to the steady boundary- 
layer flow, the steady flow must be calculated accurately. In particular, to obtain the 
large-£ asymptotic eigensolutions of the LUBLE, the base flow far downstream must be 
calculated correct to 0{w~ z ). Far downstream, the asymptotic expansion of the base 
flow is given by 


#~F + 'fi+a!s|^ + a + ifi3s!f + 'a +0 


w 




w 






(£)■ 


(D. 1 ) 


where the fractional power 72 ~ 3.774 arises from the second eigensolution of Libby & 
Fox (1963). Here F[rf) is the Blasius solution and the Pi(r}) are determined by 


where 


and 


£i(Pi) = di , 

(D-2) 

A(P) = if’ + Ff+FV, 1 

C 2 {p)~C 3 {p) = p" / + Fp" + 2F'p , -F"p, } 

C 4 (p) = £ 5 (p) s jT + Ff + 3FV - 2F"p, J 

(D.3) 


di = 1 — F' 2 , 

d 2 = 0, 

d 3 = [F /2 - 1 - 2(p2F" - p' 2 F')} - m 2 [F' 2 - 1 + 2 F'p\ + tf], 
di = P 1 P 2 ~ 3 ^ - 2 F%, 

d s = P1P3 - 3 p\p' 3 - 2F% + 2(F'p’ 4 - F"p 4 ) + 2(p , 1 p , 2 - p"p 2 ) + 2F% - 

3(1 - F' 2 ) - p 2 (F' 2 - 1 + 2F'p\ +pf). 


These equations must be solved subject to the boundary conditions Pi(0) = p^O) = 0, 
and the matching condition p' — ► 0 exponentially as 77 — > 00. The latter condition is 
necessary to ensure that vorticity decays exponentially at the outer edge of the boundary 
layer. 

The function pi(p) is the response of the mean boundary- layer profile to the 0(w~ l ) 
term of the mean pressure gradient (4.28). An exact solution for pi(p) can be obtained in 
terms of the Blasius function F(p); since the expression is cumbersome and only the wall 
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shear stress p"( 0) appears in our final results, pi(r/) is not presented here. The second term 
in the mean pressure gradient is proportional to w~ 2 . However, if the series (D.l) for <p 
were assumed to proceed only in inverse powers of w, the solution that would be obtained 
for the 0(w~ 2 ) term would not exhibit exponential decay of vorticity at the outer edge 
of the boundary layer. Exponential decay of the vorticity is restored by including a term 
of 0(w~ 2 logic 2 ) in (D.l), as discussed by Van Dyke (1964) for the symmetric flow case 
(p = 0). The equations governing p 2 (77) are homogeneous but admit an eigensolution, 
P 2 = Ai(tjF' — F). The coefficient Ai is then determined by a secularity condition on the 
differential equation £ 3 (p 3 ) = d 3 governing the 0(w~ 2 ) term. Multiplying both sides of 
this equation by F and integrating, after some manipulation we obtain 

(1 - p?) /log (t|) - j} + S l°° F(2F'p\+p\ 2 )d n = 0, (D.5) 

where Fq = F"(0) = 0.4696 and 7 = 0.33054, a constant appearing in the large-77 
asymptotic form of F(t]). Thus, setting A\ = A\+ p 2 Ai, the numerical values of these 
constants are found to be 

A t = 0.60115, Ai =-2.1656. (D.6) 


Next consider the solution for 723(77). The operator £3 is identical to £2, so that 
P 3 (r}) also contains an eigensolution B\(r)F' — F). However, in this case the function 
w~ 2 (r)F' — F) is an eigensolution of the partial-differential perturbation equation (in 
fact, it is the first eigensolution of Libby & Fox (1963)). Therefore, the coefficient B\(p) 
remains undetermined in the large-ic analysis. In §5, By{p) is determined by comparing 
the large-ic asymptotic expression for the wall shear with results obtained by numerical 
integration of the mean flow equations starting from the stagnation point w = 0. 

The ordinary differential equation governing p 3 is inhomogeneous, so that a particular 
solution is also required. The need for numerical solution for each value of p can be 
avoided by noting the form of the dependence on p of the inhomogeneous term d 3 . 
Similar remarks apply for p 4 and p 5 . Thus we set 


p 3 = Bi(r]F' -F)+p 3 + p 2 p 3 , ] 

p 4 = Aip 4 , \ (D.7) 

Ps = Bip 4 + p 5 + P 2 P5, J 

where the functions p 3 , p 3 , p 4 , p 5 and p 5 satisfy third-order ordinary differential equations 
which are independent of p. Since £3 admits an eigensolution, we ensure uniqueness of 
p 3 and p 3 by enforcing the additional conditions p^(0) = p 3 (0) = 0. The solutions for p 3 , 
p 3 , p 4 , p 5 and p 5 are obtained by numerical integration. The values of p 4 (0), 725(0) and 
725(0) are determined by a shooting technique. 

Finally, the wall shear, Uq(w ) = 4> m { V = 0), which appears in (4.29) has the asymp- 
totic expansion 


KM .F-U + frA*#.* 


H 9 + pA 1 jo 


- logic 2 


+ 


u r 


w° 


(D.8) 


Bijo + ko + p 2 l 0 , ri , _ 72 ^ 
p h 0(w 72 ) 


w J 
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where the numerical constants are calculated to be 


io =K(°)/ f o'= 3.591 , j„ =fii(0 )/FS= 0.814, 1 

h =fl(0 )/f? = 8.230, l„ s#j[(0 )//?= -3.235./ 1 ’’ 

The expression (4.31) for the curvature of the boundary-layer profile at the wall, Uq(w ) = 
<£^(77 = 0), which also appears in (4.29), immediately follows from (D.l) and (D.4), since 
p"'(0) = di{0). 
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